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Accelerator Neutrino Beams 

Sacha E. Kopp, 
Department of Physics, University of Texas at Austin 

Abstract 

Neutrino beams at from high-energy proton accelerators have been instrumental 
discovery tools in particle physics. Neutrino beams are derived from the decays of 
charged vr and K mesons, which in turn are created from proton beams striking thick 
nuclear targets. The precise selection and manipulation of the jK beam control the 
energy spectrum and type of neutrino beam. This article describes the physics of 
particle production in a target and manipulation of the particles to derive a neutrino 
beam, as well as numerous innovations achieved at past experimental facilities. 
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1 Introduction 

Neutrino beams at accelerators have served as laboratories for greater understanding of 
the neutrino itself, but also have harnessed the neutrino as a probe to better understand 
the weak nuclear force and its unification with the electromagnetic force, the existence of 
strongly-bound quarks inside the proton and neutron, and the recent revelations that neutri- 
nos undergo quantum mechanical oscillations between flavor types, a strong indication that 
neutrinos have non-zero mass. Excellent reviews of these topics are available, for example, 

in gilMllIDSllHillSSllinilMl- 

The present article discusses so-called conventional neutrino beams, those in which a 
high-energy proton beam is impinged upon a nuclear target to derive a beam of pion and 
kaon secondaries, whose decays in turn yield a neutrino beam. Such beams have been 
operated at Brookhaven, CERN, Fermilab, KEK, Los Alamos and Serpukhov, and new 
facilities at Fermilab, J-PARC, and CERN are underway. The present article cannot be 
taken as a complete catalog of every facility. Rather, the intent is to discuss some of the basic 
physical processes in meson production in a nuclear target, the manipulation (focusing) of 
the secondary beam before its decay to neutrinos, and the measurements which can validate 
the experimentally-controlled spectrum. As such, it is useful to refer to earlier papers in 
which such ideas were first developed, in addition to "state-of-the art" papers written about 
contemporary facilities. These notes will not cover so-called "beam-dump" experiments, for 
which very thorough reviews are already available |148[ 1208) . 

There are two valuable references on conventional neutrino beams to which readers may 
refer: the first is the proceedings of three workshops held at CERN P]-|S! at a time in 
which the accelerator neutrino beam concept was in its infancy. The second is a set of work- 
shops held at KEK, Fermilab, and CERN. Initiated by Kazuhiro Tanaka of KEK, this 
workshop series arose at a time of renaissance for the neutrino beam, when "long-baseline" 
neutrino oscillation experiments required new developments in accelerator technology to de- 
liver intense neutrino beams across distances 200-900 km through the Earth. Both workshop 
series are valuable documentation of the ingenuities of the experimental teams. 



2 




Figure 1: Plan view of the first accelerator neutrino experiment. Taken from Ref. |51] . 

While much has been written about neutrino interactions and detectors, comparatively 
little has been written about the facilities to produce these beams. In as much as these notes 
collect those references which may not be commonly known, I hope they will be helpful. 

2 Accelerator Neutrino Beam Concept 

The idea of an accelerator neutrino beam was proposed independently by Schwartz [f 191^ and 
Pontecorvo [183i| . The experiment, first carried out by Lederman, Schwartz, Steinberger and 
collaborators [S3], demonstrated the existence of two neutrino fiavors.^ Figure ^ shows their 
apparatus. In brief, a proton beam strikes a thick nuclear target, producing secondaries, such 
as pions and kaons. Those secondaries leave the target, boosted in the forward direction but 
with some divergence given by production cross sections d'^a/dpTdxp [pt is the momentum 
of the secondary transverse to the proton beam axis, xp ^ Pi/Pproton ratio of the 

secondary particle's longitudinal momentum along the beam axis to the proton beam mo- 
mentum). The mesons, permitted to drift in free space, decay to neutrino tertiaries. In the 
1962 experiment, the drift space was ~ 21 m. Shielding, often referred to as the "beam stop" 
or "muon filter," removes all particles in the beam except for the neutrinos, which continue 
on to the experiment. The 1962 experiment was a "bare target" beam, meaning that the 
experiment saw the direct decays of the secondaries, which were not in any way focused prior 
to their decay. 

The decays 7r± ^ /iz/^ (BR~100%), fiiy,, (BR=63.4%), and Kl vr/iz/^ (BR=27.2%) 

make the development of muon-neutrino beams the most profitable. While some muons will 

^It's interesting to note that the first accelerator neutrino beam was sufficiently new that the authors felt 
a need to put the word "beam" in quotations Though we perhaps should still do so today (no one has 
yet focused a neutrino), the current experimental facilities have certainly evolved in 44 years. 
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decay via /i eueiy^ in the drift volume giving rise to electron neutrinos, the long muon 
lifetime makes this source more of a nuisance background than a source to be exploited. Pro- 
posals have been made to produce an enhanced beam from Kl — *■ Tieue decays (BR=38.8%) 
|1H5| \W^\ \171\ inni HO] , though these have not been realized. Comparatively few experiments 
have utilized the residual z/g contamination in their beam [16, lOQj. Most z/g beams are 
produced from beam dump experiments |2081 1148j . as are z/^ beams arising from Dg tVj. 
decays |184t 1134] . For conventional neutrino beams, the neutrino spectra may be derived 
from the t:/K meson spectra and the kinematics of meson decay in flight. Some useful 
relations for the kinematics of ir/K decay in flight are given in Appendix 1X1 

The 1962 neutrino experiment didn't actually extract a proton beam. The circulating 
protons in the BNL AGS were brought to strike an internal Be target in a 3 m straight 
section of the accelerator and those resulting secondaries at 7.5° angle with respect to the 
proton direction contributed to the neutrino flux. A deflector sent the protons to strike the 
target for 25 /isec bursts ^84j- The idea of an extracted proton beam dedicated for a neutrino 
experiment came from CERN |115| I17!J] . There are at least three important motivators for 
the switch from the internal target to fast-extracted external beams: 

• The extraction efficiency onto the internal target was not perfect (about 70%, according 
to [84J or 50% in |141j . to be compared with nearly 100% for the f ast-ext r act ion [TTK] ) . 

• The CERN team developed a lens j200j to better collect the pions leaving the tar- 
get, which was much more efficient than taking those few secondaries at 7.5° to the 
beam direction, and this lens system is large (couldn't be located in or around the 
synchrotron) . 

• The van der Meer lens is an electromagnet sourced by a pulsed current which required 
short beam pulses (< 1 msec) to avoid overheating from the pulsed current. 

The second BNL neutrino run did use an extracted beam jHSllHSl, though still no focusing 
of the secondaries |146j . Extracted beams are the norm in today's experiments. 

The short beam pulses from single-turn extractions are one of the advantages of accelera- 
tor neutrino experiments: with the beam only live for a duty fraction ~ 10~^, the experiment 
(provided it has fast enough electronics) has the ability to close its trigger acceptance around 
a small "gate" period around the accelerator pulse, reducing false triggers cosmic ray muons.^ 

Neutrino experiments require expansive numbers of protons delivered to their targets. 
The 1962 experiment received 1.6 x 10^ "pulses" at an average of 1.9 x 10^^ protons-per-pulse 
(ppp)[inni- Today's experiments require 10^° — 10^^ protons on target (POT)^. Since the 
number of pion and kaon secondaries per proton grows with the incident proton beam energy, 
a good figure of merit is (POT x Beam Energy). Figure |21 shows Joules per experiment since 
the first accelerator neutrino experiment. Forthcoming experiments such as CNGS, JPARC 
and NOz/A are not shown, but are another order of magnitude in accumulated dose. 

^The BNL experiment even reduced the 25 /isec period further by gating on the RF structure of the 
circulating beam consisting of a train of 20 nsec "buckets" of protons separated by 220 nsec |109j . 
■^A LANL experiment received an impressive 10^^, but at 800 MeV/c momentum. 
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("^ pion and kaon peaks in the momentum-selected channel tunable WBB energy spectrum. 



Table 1: Tabulation of neutrino beam lines at high energy proton synchrotrons. 
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Figure 2: Compilation of total protons-on-target times beam energy per proton delivered to 
neutrino experiments at various laboratories, by date of publication. Experiments running 
concurrently in the same neutrino line are not plotted separately. 



3 Production of Hadrons in the Target 
3.1 Introduction 

Neutrino experiments require information about the production of 7r+, tt^, K~^, K~ , and 
Kl- Further, production yields, (fN/dpdfl, as a function of the secondary's momentum and 
angle emerging from the target are necessary: the secondary's momentum is related to the 
resulting neutrino energy (see Appendix^, and the production angle relates to how well the 
secondary points along the direction of the desired neutrino beam, or to the degree to which 
the secondary is captured by the focusing system. Models of secondary production have been 
derived by fitting and interpolating experimental data on p + A tc^X ot p + A ^ KX. 
The prediction of the neutrino flux starting from the yield of secondary hadrons from 
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Figure 3: Pion secondary produced by a proton striking segmented target, with definition 
of momentum components. 
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Figure 4: Fluka calculations of (left) and (right) xp distributions of tt^ in p+C 

collisions at incident momenta po = 10,20,40,80,120,450 GeV/c and pz > 0.5 GeV/c. A 
94 cm long target, 6.4 x 15 mm^ transverse size, is assumed. Taken from |178j . 



a target is the bane of every neutrino experiment. ANL, for example, performed a "beam 
survey" of the yield of secondaries from 12.5 GeV protons on thick targets of Ai and Be |149j . 
only to be surprised [90j by their neutrino flux being off by a factor of two compared with 
subsequent but more limited beam surveys I156j . The experiment scaled up the older, 
more complete d'^N/ dpdVL results to agree with the normalizations of the newer experiment 
(such was suggested by Sanford & Wang, who had tried a fit to all invariant cross section 
data |189] ) and quoted |141j 30% errors on the neutrino flux as a result. Another round 
of beam surveys was done [H^ which fixed the normalization problem and covered the full 
phase space, and these results were used in subsequent papers [7T1 1155] . As the authors of 
[7T] put it: ''The calculation of the v flux ... requires a detailed discussion, which we will 
defer to a subsequent publication. " These are hard experiments to get right. 

Figure E] demonstrates one of the aspects of hadron production predicted by Feynman 
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Pt (GeV/c) 

Figure 5: Fluka calculations of (left) pt 
various incident proton momenta Pq; (right) p' 
collisions for various values of Xp. Taken from 
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spectra of 7r+ produced in p+C collisions at 
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scaling jin4j of relevance for neutrino flux predictions. Shown is the distribution of pz and 
of Xi? ~ Pz/po for 7r+ produced by protons striking a graphite target as estimated using 
the Fluka-2005 Monte Carlo code, where po is the primary proton beam momentum 
and Pz the longitudinal momentum of the secondary (defined in Figure E}. Distributions are 
shown for incident proton momenta po=10, 20, 40, 80, 120, 450 GeV/c. The shapes of the 
xp distributions are quite similar, indicating that the pion momenta scale with pq. It is also 
of note that the integrals of these curves, i. e. the mean number of vr"*" produced per proton 
on target, grows nearly linearly with po (see Table Ej). 

Figure El demonstrates another important aspect of hadron production: the Fermi mo- 
mentum of partons inside the nucleons being ~ hc/1 fm?«200 MeV, and the fact that momen- 
tum components transverse to the boost direction are invariant, implies that the production 
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(pt) (MeV/c) 


K/ti 


10 


0.68 


389 


0.061 


20 


1.29 


379 


0.078 


40 


2.19 


372 


0.087 


80 


3.50 


370 


0.091 


120 


4.60 


369 


0.093 


450 


10.8 


368 


0.098 



Table 2: Fluka jl06j predictions for 7r+ production above > 0.5 GeV/c in a 6.4 x 15 x 
94 mm^ graphite target per incident proton. Shown are the mean number (n^r) produced 
per incident proton, mean transverse momentum (pt) for vr"*", and ratio of K/tt yields for 
several incident proton momenta po- To good approximation, (n^) oc {po)^'"^- 



spectra in transverse momentum px should be independent of xp, i.e. 

(Pn 

f{xF)g{pT) 



dxpdpj 



and the peak transverse momentum is of order 250 MeV for the secondaries. Figure El shows 
very little evolution of the pt shape for different incident momenta pq or exiting pion momenta 
Pz. That Pt does not scale (very much) is important because the transverse momentum is 
what controls the divergence of the secondary beam: mesons with = are directed along 
the beam line, and their neutrino daughters tend to follow the secondaries' direction. It is 
fortunate that the amount of pt to remove by focusing (see Section does not grow rapidly 
with pion momentum. 

The linearly increasing secondary yield with incident beam momentum has an important 
impact on neutrino beam design. It is often argued that to produce a lower-momentum 
neutrino beam one must deliver a lower momentum proton beam at the target, the rationale 
being that at lower energy machine can be operated at higher repetition rate. However, 
a given neutrino beam energy is achieved by focusing a particular secondary beam pion 
momentum. As shown in Figure |3] and Table |2l the yield at a fixed momentum appears to 
drop (approximately linearly) with decreasing proton beam momentum. Thus, the benefit 
of a lower-momentum, higher rep-rate, accelerator is cancelled by the lower pion yield per 
proton on target. The only reason for changing the accelerator energy might be to achieve 
higher secondary momenta than accessible at a lower-energy machine. 

While the above discussion of scaling is qualitatively correct, current experimental data 
indicate that these scaling behaviours are not exact. In fact, the Fluka Monte Carlo shown 
in Figures El and El being tuned to such data, demonstrates such scaling violations. 

The geometry of the target is of particular note for prediction of the neutrino spectrum. 
The geometry's significance arises because secondary particles exiting the p+A collision have 
greater probability of reinteraction in the target material for longer pathlengths. Secondary 
interactions are expected to decrease the yield of high-energy particles and increase the 
yield of low-energy particles, as reflected in the Fluka calculation of Figure IHl Plotted are 
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Figure 6: Fluka calculations of the fraction of tertiary 7r+ production from reinteractions 
in a graphite target 6.4x15 mm^ in transverse size as a function of primary beam momentum 
Po- (left) The reinteraction fraction is plotted for a 2.0 interaction length target for vr"*" with 
Pz > 0.5,5.0, etc. GeV/c momentum threshold, (right) The reinteraction fraction is plotted 
for targets of 0.5, 1.0, and 2.0 interaction lengths. Taken from |178j . 



the fraction of the 7r+ which are not produced by the primary p+C collision, but instead 
by subsequent reinteractions of the exiting particles. As also shown in the figure, such 
reinteractions occur with greater probability in high energy proton beam experiments. For 
very high energy neutrino beams, produced from high-momentum secondaries, the target is 
segmented as shown in Figure El with ~ (1 — 10) cm "slugs" separated by gaps so as to 
permit small-angle, high-momentum secondaries to escape the target with less path length 
for reinteraction. For low-energy neutrino beams, derived from low-momentum secondaries, 
such segmentation is not advantageous from the point of pion yield. ^ At CERN's CNGS 
line, the target consists 13 slugs of 10 cm graphite separated by 9 cm, appropriate for its 

^Segmented targets are of benefit for all high-power neutrino beams, however, for reducing longitudinal 
stress accumulation in the target due to heating from the proton beam. Solid targets have failed under the 
shock wave which propagates along the target's length |214j . 
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[METERS) 

Figure 7: Layout of the Cho et al. |HI] spectrometer measurement of particle production at 
ANL. The extracted proton beam is directed at a target, and secondaries are bent toward a 
Cherenkov detector by a set of dipoles. Quadrupoles keep the secondary beam focused, and 
shts or colhmators aid in the secondary momentum definition. 

focus on collection of 40 GeV/c pions [62]. 

When lacking the hadron production data which reproduces the exact conditions in a 
neutrino experiment, experimenters must rely on models to extrapolate such data to condi- 
tions of relevance for a given accelerator neutrino beam. Some of the factors which must be 
accounted for are: 

(a) interpolating between a sparse set of measurements at fixed values of secondary mo- 
menta p or transverse momenta pxi 

(b) extrapolating from measured yields off one nuclear target material to the one of rele- 
vance (Be, A£, C, W, ...) for the neutrino beam. 

(c) extrapolating to the correct projectile momentum po on the target 

(d) extrapolating to the correct target dimensions from those used in a hadron production 
experiment. 

Neutrino experiments have tried to derive these yields in auxilliary particle production ex- 
periments, either controlling or correcting for effects (a)-(d). 

3.2 Hadron Production Experiments 

Table 01 summarizes several of the hadron production experiments conducted over a range 
of incident proton momenta from 10 GeV/c to 450 GeV/c . As can be seen from the table, 
many cover only limited ranges of Xi;' and px-, owing to the geometry of the experiment. 
There are two main types of experiments: single-arm spectrometers (shown schematically in 
Figure [7j) and full-acceptance specrometers (shown schematically in Figure |HI). 
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Figure 8: Layout of the Alt et al. (NA49) [20] full-acceptance spectrometer at CERN. 
The extracted proton beam is directed at a target and the secondaries are tracked by time- 
projection chambers embedded in analyzing magnets. 



Single-arm spectrometers direct secondary particles within a small angular acceptance 
AVt into a magnetic channel in which dipoles define a secondary momentum bite Ap and 
quadrupoles are used to focus the secondaries within this momentum bite into the analyzing 
channel. Particle identification is accomplished by either TOF or Cherenkov systems or both. 
The measurements are conducted with slow-spill beams to enable single secondary particle 
counting. Normalization uncertainties on yields range from ~ (10 — 25)% due to the difficulty 
in proton counting: current-integrating toroids function well only in fast-spill (/isec) beam 
pulses, and in slow spill beams the proton intensity must be monitored either by secondary 
emission monitors (SEMs) or by the induced radioactivation in thin foils placed upstream 
of the target. SEMs are difficult to calibrate due to the decreasing secondary-electron yield 
after prolonged exposure ll(J2j , and foil-based normalizations require knowledge of the 
production cross sections for the radionuclides, which are not typically known to better 
than 10% |3I].^ In addition, spectrometer measurements require accurate knowledge of their 
acceptance. Ratios such as /ti~ or /ti^ are often better-measured. 

^Some neutrino experiments, which are fast-extracted and so could use current toroids to measure protons- 
on-target, would take their normalization from foil activation techniques, to better match what the hadron 
production experiments did for proton normalization |47j . 
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Figure 9: The parameter a defined in Equation |21 Taken from 

Full-acceptance spectrometers are a relatively recent and quite sophisticated undertaking. 
A wide acceptace tracking device, such as a time-projection chamber (TPC) is placed down- 
stream or even surrounding the target. Analyzing magnets surround the tracking system. 
For small-angle particles, downstream drift chamber planes are used. Particle identification 
is achieved by dE/dx in the tracking chamber or by downstream TOF or Cherenkov coun- 
ters. The first attempt at such a full-acceptance measurement was at CERN, in which a 
replica Cu target for the CERN-PS neutrino line was placed inside the Ecole Polytechnique 
heavy liquid bubble chamber jl72j . More recent examples are the NA49, HARP, and MIPP 
experiments, all based on TPCs. 

K'^ production is important for accurate calculation of the Ve flux from Kl neVe 
decays. While not focused, the Kl do contaminate most beam lines. ^ While Kl production 
can be approximated as 

N{Ks) = N{Kl) = ^{Nk+ + 3- Nk-) (1) 

from quark-counting arguments |59j. direct data for comparison is limited to [^Hl and |193j . 

Extrapolation must sometimes be done from a dataset collected on one nuclear target 
material to the target material relevant for a neutrino experiment. Data on pp collisions at 

^The NuTeV experiment explicitly tried to reduce this background by targeting their proton beam at an 
angle with respect to the beam line. Dipolc magnets swept the desired tt'^ and secondaries toward the 
decay tunnel, leaving the Kl to travel in the forward direction off the target |56|. 
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Figure 10: Measurement of the tt"'" yield from Be targets as a function of the target length. 
Taken from 



Po = 19.2 GeV/c [III, Po = 100 GeV/c and po = 158 GeV/c ^ are quite complete in 
xp and pt and are relevant for this purpose. Additionally, studies of the A dependence of 
cross sections at 100 GeV/c |1H] and 25 GeV/c [9Q, were used to show a scaling behaviour 

= aoA'^, (2) 

where a is graphed in Figure El This scaling proscribes how to extrapolate data taken at 
one target material to another relevant for a particular neutrino experiment. 

Neutrino targets are 1-2 nuclear interaction lengths so as to increase the fraction of the 
proton beam reacting in the target, hence the yield of secondaries. Many particle production 
experiments, however, by measuring invariant cross sections, must perform their experiments 
on thin (1-5)% interaction length targets. In so doing such experiments do not have any 
sensitivity to the effect of reinteractions of particles produced in the primary p + A as these 
secondary particles traverse the target. A measurement of particle production in "thick 
targets" is shown in Figure ^| The data are compared to a "naive absorption model" [SSI 

f{e,L)= exp{-z/Xp)exp{-tres{z)/Xs)i- (3) 

JO Ap 

where 6 is the production angle, L the target length, z the longitudinal position along the 
target, tres(^) the residual target thickness to be crossed by the secondary particle to escape 
the target. The three terms in the integral represent the probability that the proton does 
not interact up to z, the secondary is not reabsorbed, and the primary proton does interact 
between z and z + dz [22] ■ The data show excess particle production over such a naive model. 

New measurements with full-acceptance spectrometers are forthcoming from BNL E910 
|150| I132j which took thin-target data, from HARP at CERN which studied a replica of 
the MiniBooNE Be target jl90j and of the K2K Al target, and from Fermilab E907 which 
studied 120 GeV/c protons incident on a replica of the NuMI target |160j . 
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Keierence 


(GeV/c) 


Beam 


Target 
Material 


'^/^int 
(m %) 


Secondary Coverage 


HARP [78^ 


12 


PS 


Al 


5 


0.75 < j9 < 6.50 GeV/c, 30 < ^ < 210 mrad 


Asbury 25J 


12.5 


ANL 


Be 


4.9, 12.3 


p=3, 4, 5, ^ = 12°, 15° 


Cho [81J 


12.4 


ANL 


Be 


4.9, 12.3 


2 < p < 6 GeV/c, 0° < ^ < 12° 


Lundy[149j" 


12.4 


ANL 


Be 


25,50,100 


1 < p < 12 GeV/c, 2° < ^ < 16° 


Mariner 156 


12.3 


ANL 


Be, Cu 


10 


V =0.5, 0.8, 1.0 GeV/c, Q = 0°, 5°, 10° 


Abbot HI] 


14.6 


AGS 


Be, Al, Cu, Au 


1.0-2.0 


< p < 8 GeV/c, e = 5°, 14°, 24°, 34°, 44° 


Allaby [HI 


19.2 


PS 


Be, Al, Cu, 
Pb, B4C 


1-2 


p =6, 7, 10, 12, 14 GeV/c, 
9 =12.5, 20, 30, 40, 50, 60, 70 mrad 


Dekkers [HHl'' 


18.8, 23.1 


PS 


Be, Pb 


"thin" 


1 < p < 12 GeV/c, ^ = 0, 100 mrad 


Eichten ^ 


24 


PS 


Be, Al, Cu, 
Pb, B4C 


1-2 


4 < p < 18 GeV/c, 
17 <9 < 127 mrad 


Baker ^ 


10,20,30 


AGS 


Be, Al 


?? 


1 < p < 17 GeV/c, e = 4.75°, 9°, 13°, 20° 


Barton HE] 


100 


FNAL 


C,Al,Cu,Ag,Pb 


1.6-5.6 


0.3 <xf< 0.88, 0.18 <pt< 0.5 GeV/c 


NA49 [21] 


158 


SPS 


C 


1.5 


0.05 <pt< 1.8 GeV/c, -0.1 < < 0.5^= 


Aubert [SD! 


300 


FNAL 


Al 


76 


6 = 0.8 mrad, xp =0.083, 0.17, 0.25, 
0.33, 0.42, 0.5, 0.58, 0.67, 0.0.75 


Baker 


200, 300 


FNAL 


Be 


50 


6 = 3 mrad'', 60 < p < 370 GeV/c 


Baker ^1 


400 


FNAL 


Be 


75 


6 = 3.6 mrad, 23 < p < 197 GeV/c 


Atherton"^ 


400 


SPS 


Be 


10,25,75,125 


xf =0.15, 0.30, 0.50, 0.75, pr =0, 0.3, 0.5 GeV/c 


NA56/SPY 221 


450 


SPS 


Be 


25,50,75 


xf =0.016, 0.022, 0.033, 0.044, 0.067, 0.089, 0.15, 0.30, 
Pt=0, 75, 150, 225, 375, 450, 600 MeV/c 


Possible normalization discrepency with 156, 25, 81j. 



Possible normalization discrepency with Allaby et a/.jl7j. 
Full-acceptance spectrometer. 

They report angular variation between 2-3.5 mrad, consistent with then-running FNAL neutrino experiments. 



Table 3: Tabulation of published p + A hadron production experiments. 



3.3 Some Parameterizations and Models 



Without going into a complete list of all models, here are mentioned some models which 
have been employed in neutrino flux calculations. We shall not discuss some older mod- 
els/parameterizations, such as the Von Dardel |2()(ij used at CERN, Stefanski- White |194j 
used at Fermilab, or thermodynamic models |117j used at CERN, Fermilab. 

The merit of the shower cascade models is that they (claim to) contain all the necessary 
physics. They tend to be "black boxes," however, in that one cannot modify them to suit 
one's neutrino data. Such is the merit of parametric models. Comparing models to one's 
neutrino data is an unsatisfying way to evaluate systematic uncertainties, and in Sections U\ 
and IHl other techniques are discussed to adapt one's models to neutrino data. 



3.3.1 Shower Cascade Models 

Shower cascade models offer physics-driven descriptions of the cascade of particles initiated 
by a proton interaction in a nuclear target. These codes allow the user to describe a complex 
geometry of a nuclear target, impinge a beam into the target, and follow the progeny of the 
interactions through the target, allowing them to subsequently escape the target, or further 
scatter /interact to produce other particles. Such models therefore are critical to extrapolat- 
ing data with respect to the beam momentum, target material A, and understanding thick 
target effects. The state-of-the-art models include MARS-v.15 |157j . Fluka-2005^106j, and 
DPMJET-III jSU. Other models, such as GHEISHA [112], GCALOR [HI], Geant/Fluka 
[114j , or Geant4 appear to have discrepencies with published hadron production data in 
certain kinematic regimes |129j . MARS-v.15 and Fluka-2005 have been tuned to accomodate 
the SPY data, but not measurements from HARP, BNL-E910, or NA49. 



3.3.2 Parametric Models 

Malensek 

Malensek |154j parameterized the Atherton et al [201 data and included an extrapolation 
for different beam energies: 

/I ,^(l + 5e-^-'^) 
= Kp(l - xf) T n , J. (4) 

with separate parameter sets A,B,m,D for iT,K,p. Scaling to target lengths other than 
1.25Xint is done by the naive absorption model, and scaling to different nuclear targets is 
done using the data of Eichten et al jSHI- The formula maintains scaling, and the p^^ at large 
Pt was suggested by experimental data [21]. This formula fails to replicate the evolution of 
(Pt) with xf found by NA56/SPY 
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BMPT 



The authors of developed a new parameterization that fit not only the Atherton [SH! 
but also the NA56/SPY [22] data, the latter of which indicated the evolution of pt with xp- 
The functional form of their parameterization for for the invariant cross section is: 

{Ex—) = A{1 - xnTil + Bxn)x-/{1 + a'{xn)pT + 6'(a;«)p|)e-'^'(^«)P- (5) 

where XFt = E/ Emax is the ratio of the particle's energy to its maximum possible energy in the 
CM. frame, and the functions a'{xji) = a/x]^ and h\xpi) = a^/x^ control the scale-breaking 
of Pt- Separate parameters were fitted to vr"*", vr", K~^, K~ data, subject to constraints on 
ratios of positives and negatives oc (1 + axu)^ which have been well-measured previously. 
For application to other nuclear targets, the scaling Equation |21 from Barton |1H1 is applied. 
An improved version of the naive absorption model was developed for thicker targets. 

Sanford- Wang 

The Sanford- Wang parameterization [HHIEII] was used by the CERN PS ^Q^, FNAL- 
MiniBooNE |190j, K2K jl5^, and BNL beams: 

^ = Ap^il - p/po) exp[-^ - F9ip - Gpo cos^ 9)] (6) 

where po is the proton momentum, p is the pion or kaon momentum, 6 is the pion or kaon 
production angle, and the parameters A-H are fitted to experimental data, with separate 
parameters derived for tt^, tt^, K'^, and K~ . Such is a thin-target parameterization. Fits 
to Cho et al [81J and Allaby et al [T7] are found to be consistent with new data from HARP 
[75] within 10%, so this model appears quite satisfactory for "low-energy" thin target data. 

Wang [^T3] also published a variation of this parameterization suitable for extrapolation 
to higher energies upon the publication of Baker et al. [S2] '■ 

(PN 

= apQXpil - xp) exp[-6x^ - dpr] (7) 

which is quite similar to the original Sanford- Wang with the omission of the last term in the 
exponential and a new scaling for beam momentum. This function fit well to [32] with a 
factor 0.37 to account for the thick target, though it conflicts with NuMI data. 

CKP 

The CKP model |H2] apparently dates back to cosmic ray work, and was used in neutrino 
beam simulation for the CERN-PS beam [202J (adapted by [206 J and by BNL 

(Pa 



Ap {po -p)e 



~{p-a){b+ce) 



dpdVt 

which was said ^202j to be in good agreement with previous data [211 IS] ■ The effective 6""^^^ 
dependence comes from cosmic ray data ^H]. This is a thin-target model. 
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4 Focusing of Wide Band Beams 



The first accelerator neutrino experiment [H3 EE] was a "bare target beam," meaning that 
the proton beam was dehvered to the target, and the meson secondaries emanating from the 
target were permitted to drift freely away from the target. The only collimation or increase 
of flux is achieved by the relativistic boost of the secondaries in the forward direction. The 
first neutrino experiment at Fermilab jSH was likewise supplied by a bare-target beam. 

Focusing of the secondaries from the target is essential for increasing the neutrino flux 
to the detectors on axis with the beam line. In pion decay, the flux of neutrinos at a given 
decay angle 6 with respect to the pion direction is (see Appendix : 



where A is the size of the detector, z is its distance from the pion decay point, and 7 is 
the pion boost factor. If no focusing is employed, the pions diverge from the target with a 
typical angle 



where a typical pt ~ 280 MeV/c off the target was assumed (see Section Oj), and p^, ^ E.„ = 
'jm.jr. This angle of the pions off the target is larger than the typical angle of neutrinos from 
pion decay, ~ I/7, so is important to correct. Perfect focusing of pions should, in this simple 
model, improve the flux of neutrinos by ~ 25. 

4.1 Horn Focusing 

Simon van der Meer developed the idea of the "magnetic horn," |2()()j a focusing device 
to collect the secondary pions and kaons from the target and directing them toward the 
downstream experiments, thereby increasing the neutrino flux.^ 

The magnetic horn consists of two axially-symmetric conductors with a current sheet run- 
ning down the inner conductor and returning on the outer conductor, as shown in Figure ITTl 
Between the conductors is produced a toroidal magnetic field whose gv x B force provides a 
restoring force for particles of one sign (tt^ or tt"), and defocuses particles of the other sign, 
thus enhancing a z/^ beam while reducing 17^ background, for example. The focusing device 
is unusual in accelerator physics in so far as the particles must traverse the lens conductors, 
causing some loss and scattering of particles. Ref |2U2j is a thorough consideration of a 
various trajectories of particles through such a lens and the angles and momenta that can 
be focused by a particular horn geometry. 

Horns must withstand magnetic forces and the thermal load from the pulsed current and 
beam energy deposition in the horn conductors. Since the early 1970's, beam intensities were 

^The name is said to be given by the similarity of the horn's geometric shape to a Swiss alpenhorn. 
Panofsky ^175^ , however, called van der Meer's device the "Horn of Plenty." The name in the U.S. almost 
stuck [111]. 




(9) 



0^ ~ Pt/p. ~ {Pt)/p = 280MeV/p^ = 2/7, 



(10) 
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Figure 11: Van der Meer's schematic diagram of the neutrino horn, taken from |2U0j . (left) 
An optical source placed inside a reflective cone will result in exiting rays closer to the cone 
axis after several reflections, (right) A negatively charged particle emanating from a source 
is deflected by a toroidal magnetic field. 



high enough that these components become quite radioactive following extended running. 
Systems for remotely-handling any failed components are necessary ^Hl I1U31 l21Uj . Designs 
of horns now are quite refined and employ full analyses of the vibrations and strains on the 
horn (see, e.g. [SZl El] and the many contributions to [5]-[nj). Further, current- delivery 
systems have gone from large coaxial cables over to metallic transmission lines jl68| 1^ able 
to better withstand intense radiation fields and magnetic forces. The following sections 
consider various geometries of horns and their focusing properties. 



4.1.1 Conical Horns 

Van der Meer's original horn was a conical surface for the inner conductor |20m I2U2]. Such a 
device, indicated in Figure ^1 does a good job at focusing all momenta for a given angle of 
pion into the horn, 6^^ = r/i. To see this, note that the magnetic field of the device varies 
inversely with radius, B = fiQl/2'irr, and the angular deflection of the pion in the magnetic 
field (the "p^ kick"), in the "thin- lens approximation," is: 

Ae = — = 

p 27rr p 

where / is the horn current, p is the pion momentum, and x is the pathlength of the pion 
through the horn magnetic field region (see Figure IT^. 

Recalling that the incident pion angle and momentum are inversely related (c.f. Equa- 
tion [TUI), we have that the average incident angle for pions into the horn is Oj^j^ ?a {pt)/p- 
A focused pion is one in which ^out = 0, or in other words the pt kick cancels the incident 
angle of the pion into the horn. One sets this px kick to the average incident angle: 

AO = 
2iT pr p 

X = {Pt) — jr (11) 
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Figure 12: Trajectory of a and a vr" through a single conical horn focusing element. 
An axially-symmetric current sheet down the inner conducting shell produces a toroidal 
magnetic field between the inner and outer conductors of the horn, providing a restoring 
force for one sign of particles. 

This says that the pathlength in the horn should grow linearly with the radius of entrance 
into the horn, in other words a cone-shaped horn. The momentum cancels out of the final 
equation, implying this is a broad-band beam. 

Equation ^2 is derived in the limit of large source distance compared to the horn size, 
£ >> X, and the small angle approximation for the pion angle ^j^. For many horns, these 
approximations are not valid. Further, not all pions have the average {px) (not all incident 
angles are at the average or most likely angle 

4.1.2 Parabolic Horns 

It was apparently Budker who first conceived of a magnetic horn with parabolic-shaped inner 
conductors [57. Such a device focuses a given momentum for all possible angles of entry 
into the horn. It appears that such was conceived in 1961, and first attempted by his group 
at Novosibirsk to improve the collection of positrons from a target for an e^e~ collider [87j. 
The parabolic lens was studied for its efficiency in collecting mesons for a neutrino beam by 
a Serpukhov group [HS]) and first implemented in a neutrino beam at the IHEP accelerator 

A parabolic horn, like that shown in Figure is one whose inner conductor follows a 
curve z = ar^, with the parabolic parameter a in cm~^. The px kick of any horn results in 
a change in angle of 

p 27rr p 

where x = 2ar^ is the pathlength through the horn (for a parabolic conductor on either side 
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Figure 13: Trajectory of a tt"*" and a 7r~ through a single horn focusing element. An axially- 
symmetric current sheet down the inner conducting shell produces a toroidal magnetic field 
between the inner and outer conductors of the horn, providing a restoring force for one sign 
of particles. 



of the neck). Setting A^^ = ^^out ~ ^in ~ ^out ^ point source located a distance £ — f 

(focal length) upstream of the target is focused like a lens if ^out ^ 0' 

/ = ^P- (12) 

There arc two differences with the conical horn: (1) the parabolic horn works for all angles 
(within the limit of the small angle approximation), not just the "most likely angle" 6^-^ = 
{Pt)/p, and (2) a single parabolic horn has a strong chromatic dependence (its focal length 
depends directly on particle momentum p) . 

For the parabohc horn, the Coulomb scattering of particles through the horn conductors 
does not degrade the focusing quality for any pion momentum: considering a parallel beam 
incident on the horn, the spot size, S, at the focal point of the horn will be due to Coulomb 
scattering in the horn material: 

s = fez 

where 

13.6 mrad 

Oz — 

P 

is the typical scattering angle in the horn conductor, t the conductor thickness, and Xq the 
conductor material radation length. Thus 

S(x 

Thus, the quality of the focus is independent of the momentum, and improves with larger 
horn current, thinner conductors, lighter- weight materials with longer radiation lengths Xq, 
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or longer horns with larger parameter a. The fact that the focusing quality is independent of 
p means one can almost calculate a spectrum with simple ray tracing and require no Monte 
Carlo calculation |ll(Jj . To compensate the fact that particles entering the horn at larger 
radii traverse greater thickness of material t ^ 6yl + 4a?r^, horns are often designed with 
tapered conductor thicknesses, the neck region being the thickest.^ 

Figure IT^ is a demonstration performed by the Serpukhov group [87 of the momentum- 
focusing properties of the parabolic horn. A 130 MeV/c electron beam is injected into the 
parabolic horn off- axis from the left. After passing through the horn, the focusing causes a 
convergence of the electron rays at a distance from the horn equal to the focal length, after 
which the electron beam enlarges in size again. The beam size before the horn and at several 
locations after the horn is measured using photographic film. The circular spot indicates no 
aberrations despite off-axis injection and the measured focal length agreed with predictions. 



4.1.3 Ellipsoidal Lenses 

The authors of [87j (from Budker's Novosibirsk group) show, in addition to the proposed 
parabolic surface, a slightly-less tapered inner conductor shape which they term the "aber- 
rationless" surface. The nature of such an alternative inner conductor shape is better- 
elucidated in Ref. in which is shown that an ellipsoidal inner conductor surface is a 
better focusing device across wider angles of entrance to the horn. Such also appears to have 
been understood by Budker 

The ellipsoidal lens is again one in which the focal length / is a linear function of mo- 
mentum: 

f-'-^P. (13) 

where the / is again the horn current, and a and b are the major and minor half-axes (in 
cm) of the ellipsoid. 

As noted in the parabolic lens is derived in the "thin lens" approximation, and 
further requires a small-angle approximation for the particles' incident angles into the horns: 

^in ^o^in (parabolic approx.) (14) 



^This greater neck thickness is also beneficial for its greater strength. The neck is the location of greatest 
mechanical strain due to the magnetic force of the pulsed current. 
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Figure 15: Comparison of horn inner-conductor shapes required to focus particles assuming 
/ = 5m, //p=60 kA/GeV. The three curves show the required horn shape for perfect 
focusing and for the thin lens approximations of parabolic or ellipsoidal horns. The particle 
source is located at (2;, r) = (0, 0). The horns compared are "half-lenses" with vertical current 
sheet at z = 5 m. 

where = = p(GeV/c)/(6 x 10-^/(kA)). Given that ~ 200 - 400, such is more 
restrictive than the small-angle approximation required for the ellipsoidal lens: 

« 1, and ^o^in << 1, (ellipsoidal approx.) (15) 

so that the ellipsoidal lens achieves an exact momentum-focus across a wider angular spread. 
As can be seen in Figure UHl the parabolic lens is an approximation of the ellipsoid surface 
for small-angle particles. 

4.1.4 Magnetic Fingers 

Palmer |173j proposed a variant of the magnetic horn which he dubbed "magnetic fingers." 
His variation required an axially symmetric pair of pulsed conductors, but considered inner 
conductor shapes other than conical surfaces. Following numerical calculations, his inner 
conductor shape reminded him of a human digit, shown in Figure ^1 Such shapes were 
adopted for two-horn beams at BNL [ZSHZHI, and the BNL horns subsequently informed the 
designs for KEK ^T^, MiniBooNE|Tin], and JPARC [T^ . 

The numerical calculation of ideal focusing for a particle of momentum p is detailed 
thoroughly in [HHI, and dispenses with both the small-angle and thin- lens approximations, 
computing the curvature of a particle through the lens itself to obtain the required incident 
coordinates {z, r) at which a particle of momentum p should enter the horn in order to be 
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Figure 16: Inner conductor of the ideal-focusing "magnetic finger" of Palmer |173j . 

focused (^out = at the horn exit). Figure ITHl shows such a horn shape in comparison to the 
ellipsoid and parabolic approximations, in which it is assumed that the horn is a "half-lens," 
i.e.: one in which the conductor is tapered upstream of its neck, but whose current sheet 
becomes vertical at z = f, where / is the focal length (similar assumption to Figure IT^. 

The ideal surface in Figure HSl has the visual appearance of a lopsided ellipse, similar to 
lenses described by Budker's group fo^ . 

4.2 Multi-horn Systems 

Palmer |173j noted that multiple focusing elements can improve the neutrino flux because 
subsequent focusing elements can be used to "rescue" pion trajectories improperly focused 
by the first focusing element. Such a multi-lens system was adopted at CERN PS neutrino 
beam 12 7| \m\ and nearly every WBB since (see Table P). A double horn system was 
also implemented for the CERN Antiproton Accumulator |2U3j . 

Palmer |173j gives a clear motivation for the multiple lenses: a lens provides a definite 
"Pt- kick" given by A6 whose value can be calculated given the horn shape, current, and 
the particle momentum p. The horn is tuned to give a pt kick equal to this most probable 
entrance angle = {pt)/p into the horn: 

Ae = ^in. 

Many particles emerging from the target will have a pt not equal to the mean (pt), resulting 
in particles, at the same momentum p, entering the horn at a variety of angles. Assume we 
would like to focus all particles between 9^^^ = and = 29^^^. A particle entering the 
horn at will thus emerge from the horn with outgoing angle 

^out = IV - ^^1- 

A particle entering the horn with 9^^ = 9^^^ will exit at 9^^^ = 0, while a particle entering 
the horn at either ?s or 9^-^ = 26' will emerge with an angle ^q^^^ = ^j^. A particle 
beam entering the horn with angular divergence 2^j„ will emerge with divergence 9\„. 
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Figure 17: Two-lens focusing system: a second lens, significantly further from the target 
than the first, improves the collection efficiency of particles over-or underfocused by the first 
lens. The horns shown are for the Fermilab NuMI line JU] • The scale transverse to the beam 
axis is 4x the scale along the beam axis. 

A second lens far from the first will see a point source of particles with a span of angles 
to 9^^. It would be likewise expected to halve the divergence of the beam. Its inner 
aperture should be larger so as to leave unperturbed those particles already well-focused by 
the first lens. A third lens could similarly be expected to bring the overall divergence down 
a factor of 8, but must be located even further downstream to continue the point source 
approximation for the incoming particles. Techniques for design of multiple lens systems, 
including lens sizes, focal lengths, and inter-lens distances, based upon transfer matrices 
have been developed in PH]. 

A three-lens system was adopted for the 1967 CERN run|2ni I2Z1 1177j . with the second 
horn 15 m from the horn-1 and the third ~ 35 m from the target, more than half-way down 
the 60 m decay path (see Figure ITH|) . 

Serpukhov adopted a three- horn beam jSHlEHl; which had the distinction of a two-lens 
horn, shown in Figure E| the first horn consisted of two tapered regions with two "necks," 
giving the equivalent of a pair of lenses. In this sense the IHEP beam was actually a four-lens 




Shield 

Figure 18: Three-horn focusing system employed at CERN for the 1967 neutrino run, taken 
from |2()4j . Each successive downstream horn is larger to capture errant particle trajectories, 
and each has a larger inner aperture to leave un-perturbed those particles well-focused by 
the upstream horns. 
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Figure 19: Double-lens horn from the IHEP beam, taken from 1. -inner conductor, 2.- 
neck, 3.-flange, 4. -insulating ring, 5.-clamp, 6.-flange, 7.-half ring, 8. -transmission stripline, 
9. -current distribution ring, lO.-outer conductor, 11. -insulation, 12. -air cooling slot. 

system (see Figure 1301) Horns for the more recent beam lines are shown in Figure 1201 

Figure |^ shows the predicted neutrino spectrum from the two-horn system of NuMI 
at FNAL. Also shown are the components of this spectrum corresponding to the different 
pion trajectories of Figure [T7| As the angle of the neutrino parent decreases, one expects its 
momentum p ^ {pt)/^ to increase. The pions focused by only horn 1 give softer neutrinos 
than those focused only by horn 2. It is of note that the peak of the neutrino energy spectrum 
comes from particles which pass through the focusing system, while the "high energy tail" 
comes from particles which pass through the field-free apertures of the horns. Figure |221 
shows the two components from vr and K decays common to horn-focused beams. 

The NuMI beam at Fermilab implemented a "continuously variable" neutrino energy 
capability by mounting the target on a rail drive system that permits up to 2.5 m travel 
along the beam direction |139j . The target's remote control permits change of the neutrino 
energy without unstacking of the shielding elements. The utility of such a system is that it 
can assist in understanding detailed systematics of the neutrino energy spectrum observed 
in the detectors |161j . The principle of the variable energy beam relies upon Equation 1121 
since / oc p, the momentum at which point-to-parallel focusing is achieved will increase 
as the source distance is increased. Thus in the thin-lens approximation one expects linear 
dependence of the peak focused neutrino energy upon the target position Such is borne out 
by simple Monte Carlo calculation (see Figure 1221), ^"^^ by observation in the NuMI/MINOS 
neutrino data |161j .^° Further discussion is found in Section |H1 



double- neck conical horn was attempted at Fermilab. |168j but this was replaced in favor of a single 
horn [116 . 

^°The authors of Ref. and 07] note that variations of the horn current and the target positions can be 
used to vary the neutrino energy. However, these groups appear to have employed only variations in current, 
and with a goal of increasing neutrino event rate at the detectors. The authors of [931 IIU) note that linear 
lenses permit different target-horn placements to obtain low-, medium-, or high-pass beams. 
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Figure 20: Scale drawings of magnetic horns from the KEK, FNAL-NuMI, FNAL- 
MiniBooNE, and CERN-CNGS neutrino beams. All but FNAL-MiniBooNE were multi-horn 
beam lines. Adapted from |18tij . 
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Figure 21: Neutrino spectrum from the two-horn beam at the NuMI facihty at FNAL. 
The components of the spectrum correspond to the different possible pion trajectories of 
Figure El Taken from |178j . 
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Figure 22: Neutrino spectrum from the two-horn beam at the NuMI facihty at FNAL, 
showing separate contributions from tt"*", , and /i~ decays. The three graphs are from 
three settings of the NuMI beam line designed to give a low (left), medium (middle), or high 
(right) energy tune (see text and Figure EHj) . Taken from |178j . 
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Figure 23: (left) Neutrino energy spectrum at the detector hall 1040 m from the NuMI target 
for several positions of the target upstream of the first horn. Each spectrum consists of a 
focusing peak (fitted with a gaussian curve) and a high-energy tail from unfocused parents, 
(right) Peak neutrino energy (from the fitted gaussian) as a function of the target location. 
Taken from |178j . 

4.3 Quadrupole-Focused Beams 

Quadrupole-focused beams are generally less efficient than horn focusing, but they are rel- 
atively inexpensive and simpler to design, relying on magnets for conventional accelerator 
rings and they need not be pulsed, permitting use in slow-spill beams. 

4.3.1 Quadrupole Triplet 

While a single quadrupole magnet acts like a focusing lens in one plane and a defocusing lens 
in the other, pairs of quadrupoles act like a net focusing lens in both planes. Quadrupole 
triplets, furthermore, help make the containment more similar in both planes [126 , 188, .731. 
The aperture of a quadrupole is typically much smaller than for a horn, but for high energy 
neutrino beams such is not a limitation: recalling that secondaries off the target emerge 
(c./. Equation irUj) with angular spread 6 = (0.300 GeV/c)/p, a quadrupole's acceptance is 
well-matched to high-momentum secondaries. Figure |^ compares the neutrino flux from a 
horn-focused and quad triplet beam at a 500 GeV/c accelerator, for example. In principle, a 
quadrupole system provides an exact focus for a particular momentum (p) of the secondary 
beam, thus the double-peak structure in Figure results from the decays of vr's and K^s 
travelling in the secondary beam at the focused momentum (p). 

Despite providing an exact focus for particles at the design momentum (p), a quadrupole 
system is actually wide-band for detectors not too far away from the source [IBS] . As shown 
in Figure ESI particles over-focused or under-focused illuminate a detector of radius _R at a 
distance D from the source. The momentum limits of the quadrupole system are deflned by 
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Figure 24: Comparison of the quadrupole triplet focused beam with the two-horn focused 
beam proposed for 500 GeV proton beams at FNAL. The fluxes for "ideal focusing" and no 
focusing are also shown. The quadrupole beam becomes equally efficient as the horn beam 
at high energy, effectively "hardening" the neutrino spectrum. Adapted from j72j . 



the "cone of confusion," those rays coming from either the real or virtual image. 

The span of over- and under-focused particles by a quadrupole system is responsible 
for the wide-band focusing. An optical source located a distance u upstream of a lens of 
diameter 2a fills the detector with those rays emanating from the real and virtual points of 
focus at v±^^: 

(D — u)a , , 

for which the focal lengths f± are 

/± = 1/(1/«±1M) (17) 
The focal length of the quad triplet shown in Figure EEl is: 

where k = O.SBo/pa, and p is the particle momentum (in GeV/c), a the quadrupole aperture, 
Bq the maximum field at the pole tip (in Tesla). Equations ITHl and ITBl can be inserted into 
Equation [T7| to determine the limits p± of the focusing. With the quads set to focus a 

^^In contrast to geometric optics calculations, here ti± > is defined. 
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Figure 25: Illustration of limits of acceptance of a quadrupole lens pair for focusing neutrino 
parents toward a detector at a distance D from the source (target). The target momentum 
(p) is focused toward the detector, while the momentum limits p± are defined by the points 
at which the real and virtual images fill the detector of radius R. Taken from |188j . 

particular momentum (p), then u is defined hj u = {p)'^a'^ / 6 {0. 09 BqL^), and the momentum 
limits are given by 



where = a/uis the angular aperture of the quadrupole for incoming particles and = 
R/D is the desired angular illumination of the beam. As seen in Table 01 increasing the 
angular aperture of the triplet decreases the "depth of focus" (the momentum bite admitted 
by the quadrupoles), in analogy with geometric optics |188j . As D — > oo, the momentum 
bite also goes to zero, so quad focusing is appropriate only for "short baseline" experiments. 

The above discussion assumes neutrinos follow exactly the secondaries' direction. At 
{p) = 225 GeV/c, the neutrino angle with respect to the pion is ~ 1/7,^ = 0.6 mrad, to 
be compared with = 2 — 10 mrad and R/D = 1 mrad considered in Table 01 



Table 4: Momenta p± of particles over- and under-focused by a quadrupole triplet channel. 
The quads are set to provide point-to-parallel focusing for {p) = 225 GeV/c. The momentum 
limits are shown for a detector of radius i? = 1 m at a distance D = 1000 m for several 
possible angular apertures 9^^^ = a/u of the quadrupole channel (see Figure 1^ . 




(19) 



(mrad) p+ (GeV/c) p_ (GeV/c) 



2 183 318 

3 195 276 

4 201 260 

5 205 252 

6 208 246 
10 215 237 
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Figure 26: Trajectories through a (a) quadrupole doublet and (b) quadrupole triplet. The 
"F" and "D" designate quads with focusing and defocusing, respectively, for the particle 
trajectories in the plane indicated by the solid rays. The dashed rays indicate the particle 
trajectories in the opposite plane rotated 90° to that shown, for which the same quads are 
"D" and "F" , respectively. 



Quadrupole triplets are used in neutrino beams because of their near-identical contain- 
ment conditions in the horizontal and vertical views of the particles' trajectories. A doublet 
of two quads of length L and focal length / = {kL)~^ will not have equal focal planes in 
both views, as indicated in Figure EHfa), i.e. incident parallel rays will converge to a focal 
plane that is different in each view. The transfer matrix for a doublet is |126j 

where the upper (lower) signs in the ± terms indicate the FD (DF) planes. The location of 
the focal plane is given hj F = — aii/a2i, so the differing an = (1 ± kL"^) terms for the FD 
(DF) planes create an astigmatism. The equal focal lengths /* = —021 = 3/2k'^L^ in each 
view guarantee only equal angles exiting the doublet for incident parallel rays (or for the 
case of a particle source emanating from a neutrino target, we would view the drawing in 
reverse: the equal focal lengths guarantee only point-to-parallel focusing for equal emission 
angles off the target). For a neutrino beam, the secondary beam emerges from the quad 
doublet larger in the DF plane than the FD plane, which poses an aperture restriction in 
one view because quads are typically symmetric about the beam axis. 

The quad triplet shown in Figure EHtb), with "F" cells of length L/2 and a "D" quad of 
length L has a transfer matrix equal in both the "FDF" and "DFD" views |126j : 

Atr^plet=( _f^2\3/Q fA+OiVkL)' (21) 
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Figure 27: Comparison of the quadrupole triplet beam (left, taken from [SS]) used at the 
800 GeV proton beam from the Fermilab Tevatron with the "sign-selected quadrupole triplet" 
(right, taken from 



The fact that the term an = 1 in Atripiet is responsible for the near identical focusing in both 
planes. As noted in [73], subsequent quadrupole cells taking on adiabatically larger apertures 
and smaller focusing field strengths, serves to extend the momentum range of containment. 

4.3.2 Sign-selected Quadrupole Triplet 

A quadrupole system by itself focuses both signs of secondaries, thus in principle equal 
fluxes of z/^ and z7^ are obtained. In experiments in which pure z/^ or z7^ beams are desired, 
sign-selection of the secondaries must be done with a dipole to sweep out the wrong sign. 
The NuTeV experiment at Fermilab employed such a "sign- selected quadrupole triplet" |56] . 
In practice, the aperture limit of the dipole, plus the lack of focusing along the dipole's 
length, limits the wide-band acceptance of such a system by a small amount (NuTeV tuned 
to 225 GeV momentum selection, FWHM about 150 GeV). As can be seen in Figure 1771 the 
sign-selection significantly reduces the wrong-sign contamination. Wrong-sign elimination is 
especially important if running in V mode because of the lower anti-neutrino cross sections. 
Another important development in the NuTeV SSQT was the ability to target the proton 
beam at an off-angle with respect to the neutrino line, thus reducing contamination in 
the beam from Kl decays. The wrong-sign and z/e contaminations are significantly less than 
in a horn focused beam (for which 17^ ~ 10% because of unfocused particles throught the 
necks and z^e ~ 1% from muon and K decays.). 
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4.4 Other focusing systems 
4.4.1 Plasma Lens 

The BNL-Columbia group pjJ 7j proposed an alternative to the horn called the "plasma 
lens." Based on an idea from Panofsky jl74j . the idea is to place a cylindrical insulating 
vessel around the beam axis downstream of the target. The vessel has electrodes pulsed at 
~ 10 kV at its end and partial atmosphere N2 or Ar gas inside. A plasma discharge with 
current densities of ~ 10^ Amp/cm^ is initiated at the outer wall and spreads throughout 
the tube. The axial current thus produces a toroidal magnetic field, much like the horn^^. 
Particles of one sign only are focused. 

Some notable differences between a horn and a plasma lens: 

• The plasma lens in principle has no hole in its center (unlike the horn). In practice, 
neutrino beams supphed by proton beams with 400-4000 kW power would probably 
find this infeasible. 

• One can control the radial distribution of current density in the plasma to "tune" the 
magnetic field. 

Assuming a uniform current density j (in Amp/cm^) of radius R along the beam axis, 
then Bg = is present for r < R. A particle passing through this region has motion 

+ kr^ = 0, e = 

where the particle momentum is in eV. The solution to the particles motion is 

r = A sin kz, A = — 

k 

where Gq is the maximum entrance angle contained by the lens. Particles are focused parallel 
to the beam axis when kz = 7r/2, setting the desired length of the column to be L = 7r/2fc. 
The maximum radius is defined by the definition of A: 

R=^smkL = ^ (for kL = n/2) 
k k 

so the current required to focus a beam of particles emitted into the lens at < ©o is 

60 ■ 

For p = 3 GeV/c and Qq = 6°, we get J = 5 x 10^ Amperes. 

The authors report a x 3 increase in neutrino flux. The operational experience gained in 
the beam is not clear; others report initial technical difficulties [HK] . 

^^In fact, van der Meer seems to have known about Panofsky's idea of an axial current |199j . 
^■^This is a somewhat realistic example given that (pr) ~ 200 MeV/c for all pions. 
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Figure 28: Schematic diagram of the plasma lens by the BNL-Columbia group [107J. 



4.4.2 DC-Operated Lenses 

The horns of various neutrino beams have been operated at pulse-to-pulse cycle times of 
0.2 s (FNAL-MiniBooNE), 2 sec (CERN-PS, BNL-AGS, FNAL-NuMI), to 20 sec (FNAL- 
Tevatron), designed to operate in conjunction with the cycle time of the synchrotron source. 
Pulsed devices are not practical at a linear proton accelerator like those at Los Alamos or 
the SNS, or a rapid-cycle machine like an FFAG, whose Megawatts of beam power could 
prove advantageous for neutrino production |142j - |145j . nor are they practical for "slow-spill" 
beam experiments. Thus, DC-operated lenses are of advantage. Only brief mention shall be 
made here. 
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Figure 29: Schematic diagram of the DC-operated "magnetic spokes" lens of Ref. |135j . 
Magnetic Spokes 

The authors of [135^ note that the required / B -di to focus pions grows oc r as a function 
of the pion production angle off the target. For a cylindrical lens, whose focusing length i 
doesn't vary with r, this criterion requires B oc r. Given that B ^q{N I) / (27Tr) , the authors 
of jl35j chose a current distribution NI = {0.5)njr'^6, where the current is achieved by 
mounting conductors on n wedged-shaped "fins" (see Figure 1^^. each with opening angle 6, 
and carrying a uniform current density j down each side of the fins. With j uniform, then 
NI oc is achieved by having the number of conductors increase as r^.^^ 

To reduce pion absorbtion, the fins number only 8, each 8°. The return winding is 
achieved by the cables returning at the outer radial edge of the fins. The authors calculations 
show a net increase over a bare target beam of a factor of 4 with a 2.5 m long magnet carrying 
20 A. Results of the calculated fields and several pion trajectories in this field are shown in 
Figure 121 

Solenoid Lens 

As has been noted by many authors {e.g. jUj), a solenoid with axis of symmetry along 
the proton beam and target direction has the effect of transforming radial components of 
momemntum into azimuthal (angular) momentum. So, while it prevents the secondary 
beam from becoming larger, it does not by itself focus the secondaries toward a detector. 
The focusing comes from producing a gradient in the solenoid field. Ref. [02] shows results of 
a tapered solenoid which produces a field B{z) = Bo/{l + az), for example. As emphasized 
in |151j . the gradient provides the focusing through conservation of canonical momentum 
(^)0 = J^['"(-P(/) + f^^)] = 0. An advantage of this lens is that it is further from the direct 
path of the beam, while a disadvantage is that it focuses both signs of secondaries. The 
solenoid focuses certain pion (hence neutrino) momenta, which can be an advantage over a 
broad-band beam |151j . 

^"'The authors mistakenly state that the kick from a magnetic horn varies as cx 1/r. While it is true that 
the horn B oc 1/r, the pathlength of the particle through a parabolic horn grows as r^, giving J B ■ di oar, 
as required. 
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Figure 30: Plan view of the first di-chromatic neutrino beam at Fermilab. Dipoles 
momentum- and sign-select the secondaries, while quadrupoles provide point-to-parallel fo- 
cusing as they head into the decay tunnel. Taken from |46j . 

5 Focusing of Narrow-Band Beams 

In many experiments it is desirable to produce fewer neutrinos with more carefully-selected 
properties: for example, wide-band horn beams have large "wrong-sign" content (z/^^'s in a 
Ufj^ beam). Or, it might be desirable to select neutrinos of a given energy for study of energy- 
dependence of cross sections or neutrino oscillation phenomena at a particular neutrino 
energy. 



5.1 Dichromatic Beam 

Fermilab was the first to pioneer the so-called "di-chromatic neutrino beam" |147j . and the 
high event rates possible yielded rapid physics results UHl ESI- Such a beam, shown 
schematically in Figure EOl uses dipole magnets downstream of the target to sweep out 
quickly wrong-sign secondaries from the neutrino channel. In the first di-chromatic beam, 
two quadrupole magnets were used to provide point-to-parallel focusing of those secondaries 
of the momentum selected by the dipoles. The monochromatic secondary beam of pions and 
kaons is sent into the decay tunnel, where they decay. Following the construction of the 
SPS at CERN, a similar dichromatic beam was built there |123j . with physics results coming 
from the CDHS[Slj, and CHARM [IH] detectors. 

The term "di-chromatic" comes from the two distinct neutrino energies produced in such 
a decay channel. The decay of a pion or kaon secondary results in a neutrino of energy 

p _ (1 - irn^/m^n,K)?)E(^,K) . . 

(1 + 72^2) ' ^^^) 
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Figure 31: (left) Graph of ^ + E^Ymi CITFR calorimeter as a function of 

transverse radius of the neutrino interaction with repect to the beam axis jlHI- (right) 
Graph of the visible energy in the CITFR for the +190 GeV/c beam setting, with the vr and 
K radial bands shown separately 



where 9 is the angle between the neutrino and meson direction, and 7 = -E(,r,if)/'"^(7r,_ft')- The 
momentum of the secondary beam is fixed, but the presence of both pions and kaons lead to 
two possible values for the neutrino energy. The possibility for off-angle decays of the (vr, K) 
beam can change E,y. Figure IHTl shows this kinematical relationship in the Caltech-Fermilab 
neutrino detector located 1300 feet from the end of the decay pipe: neutrino interactions 
reconstructed in their detector at large transverse distances {i.e.: large (vr, K) decay angles) 
from the beam central axis show a smaller total energy deposition in the detector, though 
two distinct bands are observed, arising from pion and kaon decays. 

The channel downstream of the target starts producing neutrinos as soon as secondaries 
decay. Decays before the momentum- and sign-selection are achieved result in a "wide- 
band background" under the two energy peaks in Figure |^ For this reason, the proton 
beam is brought onto the target at an angle off the axis of the decay tunnel, resulting in 
such "wide-band" secondaries decaying preferentially away from the neutrino beam's axis. 
Further, momentum-defining collimators are placed along the neutrino channel to better 
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Figure 32: Two-horn system for the BNL narrow-band beam, with colhmators and beam 
plugs to stop unwanted pion trajectories. Pion ray #1 is transmitted, while pion ray #2 is 
stopped by the collimator between the two horns. Likewise, rays passing through the two 
beam plugs would be expected to be attenuated. Figure adapted from |192j . 



eliminate off-momentum secondaries from the beam. In fact, these considerations, plus the 
upgraded capabilities of running the Fermilab Main Ring at 400 GeV/c primary momentum, 
led to an upgrade of this dichromatic beam [OEl EZj with larger primary targeting angle to 
reduce the wide-band backgrounds and better momentum selection to reduce wrong-sign 
contamination. 



5.2 Horn Beam with Plug 

The wide-band horn-focused beam, referring to Figures 1^ and fTTl produces a span of neu- 
trino energies corresponding to a variety of particle trajectories through the focusing system. 
To cut off the largest range of neutrino energies, it is desirable to eliminate those particles 
which travel through the field-free "necks" of the horns. Such was attempted at CERN 
|177| l58j by placing a Tungsten block (beam "plug" ) at the end of the usual target to help 
attenuate those high energy pions which tend to leave the target at small angles {9-,^ ~ '^/i-k)- 
The collimation for a narrow-band beam was refined in a series of experiments at BNL 
j75| [7n] , in which two beam plugs and a collimator located in between the horns were used 
to attenuate all but the desired trajectories, as shown in Figure Referring to Figure IT7| 
further eliminating those particles which do not cross the beam center line between the 
two horns has the effect of cutting all but the smallest momenta, as is achieved with the 
collimator between the two horns in Figure |221 A similar proposal was made at Fermilab 



5.3 Horn Beam with Dipole 

As noted in ^U], a dipole magnet placed in between the two horns of a wide-band beam 
has the effect of achieving better momentum and sign selection. As shown in Figure ESI a 
dump for the primary beam must in this case be placed in the target hall, just like in the 
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Figure 33: NBB achieved by a momentum-selecting dipole placed in between the two horns 
of a WBB. Taken from |lOj. 

dichromatic beam, which is somewhat of a challenge for high-intensity neutrino beams. In 
practice, the aperture restriction of the dipole does attenuate some of the pion flux. 

5.4 OfF-Axis Neutrino Beam 

The idea for an off-axis neutrino beam was first proposed by BNL experiment E889 |5(Jj . 
Many of the kinematic features of off-axis pion decay were worked out in Ref. |185j . The 
estimates of the on- axis WBB flux in Section |3] made implicit use of the fact that the energies 
of neutrinos emitted along the axis of travel of the secondary pion or kaon is linearly related 
to the meson energy. The problem of achieving a particular energy NBB thus reduces to 
focusing a particular energy meson beam. 

In the limit that mesons are focused and travel parallel to the decay pipe axis, the BNL 
E889 team noted that under some circumstances nearly all mesons of any energy could 
contribute to generating the same energy of neutrino. While Equation |221 states that the 
neutrino and meson energy are in fact linearly related for on-axis decays {9 = 0°), the 
relationship is more complex for neutrinos observed to emerge at some angle with respect 
to the beam due to the denominator. Equation |221 is graphed for several particular decay 
angles in Figure 123 

Figure OH has an interesting interpretation: for on-axis decays, the neutrino energy is 
related to the meson energy. For off-axis decays, this relationship is weaker. Thus, for 
large off-axis angles, nearly any pion energy makes about the same energy of neutrino. A 
broad-band pion beam, therefore, can be used to generate a narrow-band neutrino spectrum. 
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Figure 34: Neutrino energy from pion decay as a function of pion energy, for several choices 
of decay angle between the neutrino and pion direction. Taken from ||50j. 
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Figure 35: Neutrino energy spectra at a distance of 1 km from the proposed BNL beam for 
several off-axis angles. Taken from (SUj- 
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Neutrino Energy (GeV) 

Figure 36: Calculated flux from the NuMI beam at 2; ~ 750 m and 110 mrad off-axis angle, 
corresponding to the location of the MiniBooNE detector. Two peaks, from 71 and K decays, 
can be seen. Taken from |13j . 

The BNL team proposed such a NBB spectrum for a search for z/g oscillations, 

since NC z/^ interactions of any energy can leave small energy depositions in a detector 
which mimic z/g interactions. Thus, cutting down all z/^ energies which contribute to NC 
background is of value. They proposed placing a detector a couple of degrees off the beam 
axis for their new beam line, thereby choosing the particular NBB energy to be achieved. 

Figure EHl shows, for the beam configuration and detector distance in the BNL proposal, 
the neutrino energy spectrum from pion decays at several off-axis locations. In addition 
to the lower, narrower, energy spectrum at larger off- axis angles, it may be noted that, at 
certain energies, the flux at the peak actually exceeds the flux at that same energy in the 
on-axis case. Thus, the fact that all pions contribute to approximately the same neutrino 
energy can, in part, compensate for the loss of flux at off-axis angles, from Equation El 

The proposal, not approved, has since been adopted by teams at JPARC |128j and 
Fermilab |170j . which will employ the narrow-band off-axis beam to search for 
oscillations. The first detection of neutrinos from an off-axis beam is at Fermilab, where the 
MiniBooNE detector is situated 110 mrad off-axis of the NuMI beam line at a distance of 
~ 750 m from the NuMI target. Neutrinos from NuMI have been observed in MiniBooNE 
[T3] . The off-axis angle is sufficiently large that both peaks from vr and K decays can be 
seen (see Figure permitting use of the MiniBooNE detector to derive the ti / K ratio of 
the NuMI beam. 15 

^^The ability to resolve separate pion and kaon peaks at large off-axis angles, as well as the low systematic 
uncertainties in predicting the flux of an off-axis beam were studied in 176j . 
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Figure 37: Schematic diagram of the MiniBooNE beam hne at FNAL, taken from |61j . 

6 Decay Volumes 
6.1 Decay Tube 

Decay volumes are drift spaces to permit the pions to decay. For a 5 GeV pion, 7 ^ 35, and 
•jPcTt, ~ 280 m. This sets the scale for how long the decay tube should be if just 63% of the 
pions for a 2 GeV neutrino beam are to be allowed to decay. As noted by [202 , the decay 
pipe radius is also of importance, and has to be as wide as practical for efficient low neutrino 
energy beams: in general low-energy pions are not as well focused in a horn focused beam, 
and have a divergence which will send them into the decay volume walls before decaying. 

Decay tubes are often evacuated. The same 280 m mean flight path, in air, represents 
0.9 radiation lengths {Xq = 304 m for air at STP), and 0.26 nuclear interaction lengths 
(Knt — 1080 m). Thus, a pion drifting in air at atmospheric pressure would have a ~ 26% 
chance of being absorbed by a collision, and those that are not lost will suffer multiple 
Coulomb scattering of a typical magnitude of 2.8 mrad. Such scattering angles are already 
significant compared to the ~ I/7 opening angle between the muon and neutrino in pion 
decay, which is 14 mrad for a 10 GeV/c pion decaying to a 4 GeV neutrino. Other decay 
tubes, such as KEK [T3], are filled with He gas to reduce absorption and scattering. 

Because scattering can do much to defocus the secondary beam already focused by the 
horns, particular care is given to the entrance windows to decay volumes. The needs of 
mechanical strength for the large evacuated chamber must be balanced against placing sig- 
nificant scattering material in the beam. NuMFs 2 m diameter decay pipe has a composite 
window with a 3 mm thick, 1 m diamter Aluminum center and thicker steel annulus at larger 
radius to reduce pion scattering and heating from the unreacted proton beam. The CNGS 
beam has a thin Ti window |lOlj . 

T2K |131j has a fiared decay volume which enlarges at its downstream end, as shown 
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beams-eye view 

Figure 38: Schematic of the decay volume for the T2K beamhne, taken from |131j . 

in Figure 13^ This beamhne is envisioned to support experiments at two remote locations, 
one at Super-Kamiokande and also a future "Hyper-Kamiokande" site. It is envisaged to be 
an off-axis beam (see Section 15. 4j) of about 2-3° to both these sites. The flared beam pipe 
permits tuning of the off-axis angle (hence {EJ)) as the experiments require. 

6.2 Hadron Hose 

Fermilab proposed building a focusing device along the length of the decay pipe which 
would enhance the neutrino flux and reduce systematic uncertainties in predicting the energy 
spectrum of neutrinos |127j . Based on the "beam guide" idea originally proposed by van der 
Meer j2Ulj , the device consists of a single or multiple wires travelling axially down the length 
of the decay volume which are pulsed with ~ 1 kA of current, providing a weak toroidal field, 
but long focusing length (the full particle trajectory before decay). As indicated in Figure IH^ 
such focusing draws particles diverging toward the decay volume walls back toward the beam 
center, where they can decay without absorption on the walls. 

The hadron hose can increase the neutrino event rate to experiments by 30-50% because 
pions heading toward the decay pipe walls are drawn back toward the beam centerline. 
Improved probability for pion decay can also be achieved simply by constructing a larger 
diameter decay volume, but such is quite expensive due to the extensive shielding which 
must surround the decay volume in high-power neutrino beams. Thus, the hadron hose may 
be viewed as an active decay volume, a low-cost alternative to the large-diameter passive 
decay volume. 

The hadron hose provides a second benefit which is less obvious: the spiral orbits essen- 
tially randomizes the decay angles of the pions leading to neutrinos in the detector. This 
is beneficial for two-detector neutrino experiments, because the "near" and "far" detectors 
often observe slightly different energy spectra just due to the sohd angle difference between 
the two detectors. Recalling that the neutrino energy is Ey = (0.43£'7r)/(l + 7^6'^), high 
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Figure 39: Schematic diagram of the hadron hose focusing device proposed by |127j . based 
on the "beam guide" of van der Meer [20 IJ . Secondary particles entering the decay volume 
spiral around the current-carrying wire until they decay to neutrinos. 

energy pions which decay just in front of the near detector can result in neutrinos hitting 
the near detector for a wide span of angles 9, lowering the neutrino energy as compared to 
the neutrinos reaching the far detector at 6' ~ 0. The randomization of decay angles, caused 
by the spiraling orbits in the hose field, is discussed further in Section |H1 

The focusing might naively be expected to converge all particles into the wire, causing 
large absorptive losses of pions: pions emerge from the target in the radial direction, and 
the radial restoring force causes many pion trajectories to cross the wire. However, multiple 
Coulomb scattering of the pions and kaons in the upstream horns and entrance window to 
the decay volume leads to some azimuthal component of pion momentum, causing the pions 
to enter the decay volume and execute spiral orbits around the hose wire []i)2j, as indicated 
schematically in Figure EHl Analytic expressions for particle orbits in the hadron hose field 
have been computed |187t I162j . 

Placing a high-current wire in the evacuated decay volume poses some technical chal- 
lenges, as discussed in |127j . Namely, the wire's heat induced by PR as well as energy 
deposition from beam particles must be dissipated sufficiently by blackbody radiation, the 
wire's voltage must be shown not to break down in the heavily ionized residual gas of the 
decay volume, and the long-term tension applied to the wire segments must be ensured not 
to cause plastic flow ( "creep" ) of the wire material such that a failure occurs. 

6.3 Muon Filter 

The muon filter is the part of the beam line required to range out muons upstream of the 
neutrino detector. Keeping in mind that ~ 2 MeV/(g/cm^), and recalling for steel 
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Figure 40: Plan view of the IHEP neutrino beam. Taken from [35j. In addition to the 
four-lens neutrino channel, a small side channel at 87° could be activated for delivery of 
secondary mesons from the target to the neutrino experiments for calibration purposes. 

(often used in shielding) that p ~ 8 g/cm^, ^ ~ 1-6 GeV/m for those nuisance muons. 
The first neutrino experiment in fact had to lower the AGS accelerator energy to 15 GeV to 
reduce the maximum muon energy and thereby reduce the muon "punch-through" jHl|. The 
origninal neutrino line at Fermilab, which had an earthen "berm" sufficient to stop muons up 
to 200 GeV/c, had to be reinforced with 20 m of lead and 140 m of steel shielding following 
upgrades of the accelerator complex to run at 800-900 GeV proton energy 

The location of the upstream face of the muon filter defines the maximum pion or kaon 
drift time before decay to muon and neutrino. It is expensive to construct a decay tube 
that allows most focused pions to decay. For example, the CERN PS neutrino beam, with 
80 m decay volume, would allow 25% of pions and 90% of kaons to decay, assuming that 
5 GeV particles are being focused. In the case of NuMI, with 725 m of drift space and 
~10 GeV/c focusing, these numbers are 73% and 100%, respectively. The length of the 
decay volume also impacts the level of z/g content in the beam, since much of it arises from 
TT nUf, {eUf,Ue)Uf, decays. 

The idea of a moving beam dump (see Figure will be employed by MiniBooNE to 
demonstrate their u^s come from oscillated z/^'s from pion decay, not the "instrinsic" from 
the beam of — vrez/e decays or /i decays. The moving beam dump was first employed 
by jHl] to show their neutrino candidates were from meson decay and not from interactions 
in the shielding. Figure [T] shows a lead block that was placed close to the target to stop 
presumably all but a few K decays, and indeed the neutrino rate decreased in proportion 
to expectations. CERN's original WANE beam also apparently had a moveable mid-stream 
beam stop, but this was never employed |123j . 

Low- intensity beam lines combined the proton beam dump with the muon filter I40j . 
However, for modern beam lines, a dedicated proton dump is required because of the intense 
beam power. The NuMI beam, for example, is designed for a 400 kW proton beam, of which 
70 kW heads for the beam dump, requiring water cooling, a special Aluminum core, etc. 
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Figure 41: Plan view of the CERN PS neutrino beam from, taken from [HH]- A small 
mercury-filled tube in the muon shield could be drained temporarily, exposing the down- 
stream neutrino detectors to muons from the beam for calibration. 



Accident conditions are even more problematic: the beam stop must allow for errant proton 
beam missing the target and striking the dump directly. This is an even greater concern for 
upgrades to NuMI, CNGS, and JPARC, the latter two of which will use a graphite core. 

A common problem in muon shielding is leakage, not attenuation |181j . Staggered assem- 
bly with no "line of sight" cracks is crucial to good shielding design. This uniformity impacts 
the ability of downstream muon instrumentation to make meaningful measurements of muon 
intensity and lateral profile. Such measurements, which can provide information on the neu- 
trino flux and even the energy spectrum, are distorted by cracks which let lower-energy 
muons through, as has been observed at NuMI. 

There have been a couple clever trickes to temporarily "let down" the muon shield of 
an experiment for the purposes of calibrating the neutrino detectors with particles (muons, 
pions) of known momentum. The Serpukov beam could calibrate its spark chamber and 
bubble chamber experiments [33] using a small channel in their shielding at an angle 87° to 
the primary beam axis. Shown in Figure 1^ this channel permits secondaries from the target 
to be focused in a quadrupole-dipole system and be delivered directly to the experiments. 
Such a calibration test beam was also utilized by the NuTeV experiment at Fermilab JJL9J. 
Another trick employed at the CERN PS neutrino beam in 1967 was to install a mecury-filled 
tube which penetrated the entire 20 m muon shield |177j . shown in Figure 1^ The mercury 
from this tube could temporarily be drained, exposing the heavy-liquid bubble chamber 
(HLBC) from Ecole Polytechnique to muons at the end of the decay volume. 
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7 Flux Monitoring 



7.1 Primary Beam Monitoring 

The monitoring of the primary proton beam, as far as it impacts the physics of a neutrino 
experiment, is hmited to requiring knowledge of the proton beam just upstream of the 
target. Specifically, paramters such as the total intensity of the beam striking the target 
(both integrated over the lifetime of the experiment and on a per-pulse basis, since many 
experiments suffer rate-dependent effects), the position, angle, divergence and spot size of 
the beam as it is about to strike the target. 

The proton flux delivered to the neutrino target can be measured in a variety of ways. 
Fast-extracted beams can use current toroids, and NuMI has recently demonstrated cali- 
bration of such a device to ~ (1 — 2)% over the first year of operation using precision test 
currents. In the past, many experiments would often take their normalization from foil acti- 
vation techniques, which measured the residual activity of gold |17], M |177j . or polyethylene 
jnZj foils placed in the proton beam upstream of the target. Such techniques are typically 
precise to (5 — 10)%, due to imprecise knowledge of production cross-sections for these ra- 
dionuclides. One motivation for using such foil techniques was to better match what the 
hadron production experiments did for proton normalization jl?], but this is becoming less 
imporatant as experiments are relying on more than one hadron production experiment. 

The proton beam profile has in various lines been measured by segmented ionization 
chambers iOl 113 EEl ED] , Aluminum SEMs HHg, W wire SEMs [M!, ZnS screens|TT71 R7I 
ISEi- Many of these techniques no longer work in high-power neutrino lines: the large proton 
fiuences motivate the need to reduce beam scattering and loss along transport line, as these 
cause irradiation and damage to transport line magnets. Further, the proton beam's power 
can significantly degrade the performance of an interceptive device in the beam. At NuMI, 
the profile is measured at the target with a segmented foil Secondary Emission Monitor 
(SEM) [TTTj . 

7.2 Secondary Beam Monitors 

Instrumentation placed directly in the secondary (vr, K) beam of a wide-band beam is rel- 
atively rare, since it must cope with quite high rates and can substantially affect the neu- 
trino flux. A few notable examples exist. CERN proposed placing a spectrometer and 
Cherenkov counter system downstream of their horns to measure n/K fluxes after the horn 
focusing. |18nj . Such was a "destructive measurement," from the point of view of neutrino 
running, but would have yielded an in situ analysis of hadron production and focusing. In 
Figure 1^ this spectrometer is indicated by the Cherenkov counter just below (beam left) of 
the muon filter tilting at an angle which points back to a thin specrometer magnet (curved, 
in front of the horn R2). A test of this system was conducted |182j . though backgrounds 
from 5-rays in the beam and produced in the spectrometer appear to have been difficult. 

KEK placed a Cherenkov counter, shown in Figure |32l in their secondary beam for two 
brief periods during their run jT^ I158j . This system placed a wedge-shaped sherical mirror 
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Figure 42: The K2K pion monitor. Taken from [TK] . 



at 30° to the beamhne to direct Cherenkov hght out to a PMT array several meters away 
from the beam axis. Assuming that all particles in the beam are pions, the Cherenkov ring 
sizes provide the pions' momenta while their location on the PMT array provide the pions' 
direction off the beam axis. Substantial (~ 30%) substractions were made for electromagnetic 
particles in the beam. With the (p^r, 9^^) information, a modified fit to the Sanford-Wang 
parameterization |189j is possible. To avoid detections of protons in the Cherenkov counter, 
it could measure the neutrino spectrum above 1 GeV (pions above 2.3 GeV/c), which is 
approximately the location of the maximal flux (see Figure I3T|) . 

BNL[80j and the Fermilab NuMI beam |138j placed segmented ion chamber arrays directly 
in the secondary beam as beam quality monitors. The NuMI chambers must contend with 
~ 2 X 10^ particles/cm^ /spill and are exposed to ~ 2 GRad/yr dose, necessitating moving 
away from circuit board technology as in [HO] to all-ceramic/metal design |138j . Because of 
the large fluxes of photons, electrons, positrons, and neutrons in the secondary beam, neither 
chamber was used in a flux measurement. The BNL chambers were placed midway down the 
decay volume, while the NuMI chambers were located right upstream of the beam absorber. 
In the case of the NuMI beam, the flux at the hadron monitor is dominated by unreacted 
protons passing through the target, so the device serves as a useful monitor of the proton 
beam targeting, as well as a check of the integrity of the target. The CERN WANF beam 
[2H1 1123j had split-foil SEMs downstream of the target but upstream of the horns to ensure 
beam was on target, and the CNGS beam will do likewise |101j . 

The dichromatic beam at Fermilab had an elaborate secondary beam system which was 
crucial for making flux measurements and which enabled absolute neutrino cross sections to 
be measured. The narrow, momentum-selected secondary beam permits reasonably small- 
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Figure 43: The FNAL dichromatic neutrino beam, with secondary beam instrumentation 
indicated. An ion-chamber array and an RF cavity measure absolute particle flux in the 
secondary beam, and a Cherenkov counter measures relative abundances of e, vr, K,p in the 
secondary beam. Taken from [IBj . 



diameter instruments to be inserted or removed from the secondary beam. These detectors 
included two ion chambers which measured total particle flux, an RF cavity which was 
used to corroborate the ion chamber measurement, and a Cherenkov counter which could be 
scanned in pressure to measure the relative abundance of e, vr, p in the beam (subsequently 
normalized to the total flux determined by the ion chambers). With this system in place, it 
was not, in principle, important to know the number of protons delivered to the target in 
order to estimate the neutrino flux. 

The ion chambers were carefully calibrated. Linearity with particle flux was demonstrated 
by comparison to the proton fluence on target measured by a beam current toroid. Stability 
in time was shown by comparions to the two ion chambers' relative signals. Studies were 
done to show that material upstream of the ion chambers contributed negligible signal in 
the form of 5-rays (the ion chambers were placed well-downstream of any shielding), and 
the signal response was carefully studied as a function of relative particle abundances, since 
heavy protons cause a rise in signal due to nuclear interactions in the ion chamber materials. 

The Cherenkov counter was essential to this measurement because the vr — > /ii/ and 
K ^ fiu decays contribute to different energy neutrinos in the dichromatic beam, as shown 
in Figure 1^ A plot of the relative abundance during +200 GeV/c secondary beam running 
is shown in Figure |^ Measurement of these two individual fluxes absolutely, along with the 
known momentum bite of the dichromatic channel, allows absolute flux predictions which 
can then be compared with the event rates in Figure 1^ to derive cross sections, independent 
of knowledge of protons on target. 
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Figure 44: (left) Pressure scan of the Cherenkov counter in the CITFR dichromatic secondary 
beam, showing clear peaks in proportion to the vr"*", , and p fluxes in the secondary beam. 
The plot area is normalized to the total particle flux measured by the ion chambers, (right) 
Ratio of K/ TT fluxes vs secondary beam momentum for positive and negative beams. Taken 
from 



7.3 Muon Beam Monitoring 

There are two kinds of muon systems that have been built: flux measuring systems and di- 
agnostic systems. Flux monitors attempt to use the tertiary muons to yield a measurement 
of the neutrino flux. This is a plausible idea, since muons come from same decays as the 
neutrinos. Not all beamline geometries are conducive to such flux measurements, however, 
because of either decay kinematics or because the shielding imposes limitations on the frac- 
tions of the muon flux visible to the muon detectors. First, the muon detectors must be 
placed downstream in the muon filter, imposing a lower threshold of muon momentum and 
thereby cutting off direct measurement of the lowest part of the neutrino energy spectrum. 
Second, the solid angle acceptance of such muon detectors (even in absence of intervening 
shielding), can be prohibitively small, especially at low energy where pion decays are wider 
angle and may not intercept the muon detectors at the end of the decay volume. Diagnostic 
muon monitor systems may be similar to flux monitors, but smaller fraction of the 

muon spectrum. They must, however, be available "online" during neutrino running, which 
poses different constraints in terms of simplicity of construction and radiation tolerance. 

While the first neutrino experiments at BNL and CERN did not in any way measure 
their neutrino flux, the follow-on experimental run at BNL 166] did so, and was the first 
attempt to measure a neutrino flux using the tertiary muons. Emulsions were placed in 
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Figure 45: Emulsion measurements of muon fluxes in the steel shielding in the BNL neutrino 
experiment jHE]- The curves are flts using the CKP model jH21- 



seven "probe holes" in the steel shielding at the end of the decay region, each probe hole 
containing up to four emulsions at different transverse distances to the beam axis. Since the 
experiment had no focusing of the neutrino parents, this measurement, in the limit of no 
multiple scattering in the steel, should help corroborate their neutrino flux calculation. The 
experimenters report an error (20 — 30)% from these measurements, due to the incomplete 
phase space sampled by the emulsions in the steel and the inability to go below a certain 
threshold (steel thickness) for fear of backgrounds from upstream hadrons. Their flux data 
is shown in Figure |^ This tuning was done using a short run with a 3" Be target to allow 
comparison with the existing thin-target hadroproduction data and with a 12" Be target 
which was used during the neutrino run. 

CERN subsequently measured its neutrino spectrum using muon system measurements 
|58j . There are multiple challenges: (1) they re-parameterized Sanford-Wang from this flt, so 
have to hope thick target effects simply scale the flux, not modify it; (2) focusing effects have 
minimal effect on the spectrum; (3) muons from K decays only easily distinguished at large 
lateral offsets from beam axis, so have to assume tt /K ratio from external 'beam survey' data 
^^; (4) the flrst few meters of shielding have large hadron shower content, so one must only 

similar measurement of the neutrino flux was attempted using the muon system at the IHEP beam 
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Figure 46: Demonstration of the muon flux- fitting procedure performed at CERN to deter- 
mine the neutrino spectrum, (left) muon fluxes at several lateral locations transverse to the 
beam axis and at several longitudinal depths in the muon filter, (right) Sanford-Wang fU^H] 
parameterization of pion yields from the target (PN/dpdQ before and after the muon fit, as- 
suming all target particles are created in primary interactions and the nominal Sanford-Wang 
K/tt ratio, (right) Neutrino flux at the bubble chamber before and after the fit. Figures 
taken from (207j . 

measure the flux above some threshold or perform a significant subtraction (CERN claims 
50% for the first data point at 1.7 m, 12% in the second, and 6% in the third - these numbers 
were confirmed by measurements subsequently made with a W beam plug after the target, 
which demonstrated drops in roughly these proportions); (5) Both yU+ (proportional to z/^ 
flux) and fi~ (proportional to the z7^ flux, which CERN claims is only 0.2% of their flux); 
(6) the measured flux has to be corrected by ~6% for 5-rays (the correction was obtained 
using emulsions placed on the chambers. A 'beam survey' of Ref. ^Tj became available after 
this work, and agreed in neutrino flux prediction to within 10-15%, which was their stated 
uncertainty.^'^ 

The large flux of 6 rays produced in the muon filter shielding makes flux measurements 
in a wide band beam quite challenging jl2H IT^ . Charge- integrating detectors like solid 
state detectors or ionization chambers will measure a combination of muon signal and 6- 
ray background. As shown in the left plot of Figure El the electron component can be a 

|23| . and at the higher energies there the K/tt ratio uncertainty seems to be a bigger effect. 

""^^Unfortunately for present-day beams, this method may not be easily apphcable at high-intensity neutrino 
beams. Higher intensities have required thicker beam stops to absorb hadronic showers from the remnant 
proton beam, and the neutrino energies (hence muon momenta) of interest have decreased for long-basehne 
neutrino osciUation searches. 
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Figure 47: (left) Monte Carlo calculation of the transverse distribution of muons and knock- 
on electrons in one of the muon pits at the CERN WANF beam, (right) Measured track 
angle in an emulsion sample placed in a 100 GeV/c muon beam, showing tracks at zero 
degrees (beam muons) and wide-angle tracks (knock-on electrons). Taken from jl2j . 



significant fraction of the signal and has a different lateral shape than the muon beam due 
to multiple scattering of the electrons. CERN employed a series of emulsions to count the 
tracks as well as the scattering angle of each track (straight-through tracks are presumably 
from beam muons while wide-angle tracks are from knock-on electrons), thereby obtaining 
a /x/e "correction factor" by which future muon measurements are adjusted. Thus, absolute 
flux measurements are limited to 3-5% due to the counting in emulsions and the ability to 
separate electrons from muons. 

The technologies employed for muon monitors vary, though they must become increas- 
ingly simple and radiation-hard as the intensities increase and the ability to access these 
remote muon pits is reduced. Ionization chambers were used in horn-focused wide-band 
beams such as the CERN PS neutrino beam 5HI I177j . the BNL neutrino beam [5UI|.the 
IHEP- Serpukhov beam [^3], the KEK neutrino beam |124[ [T3]. the Fermilab NuTeV beam 
PTH] . the Fermilab NuMI beam [HHl, and the CERN CNGS |M|- Solid state muon detec- 
tors were used in the CERN West Area Neutrino Facility (WANF) beam line off the SPS 
|7^ll22tri23j . Maximum muon rates ranged from 5 x 10^/cm^/spill ^ to 5 x 10^/cm^/spill 
|138j . ANL jHO], CERN |58j, and Serpukhov had plastic scintillators in the downstream 
portion muon shield, where particle fluxes are lower ~l/cm^ /spill). 

Muon monitor systems have done real diagnostic work: the CERN muon chambers at 
the PS neutrino line detected a flux asymmetry which was eventually traced to a magnetic 
field asymmetry in the horn [95 . Also at CERN, a misalignment of some of the target 



54 



J ID 
0.1 



Monte Carlo 



J 



No backgrounds, only "clean" muons. 
No momentum smearine. 




muon momentum at 1° (GeV/c) 

Figure 48: Muon momenta seen 7° off-axis from the MiniBooNE decay pipe from tt and K 
decays. Taken from 



liall components was detected at the WANF line: they could achieve higher muon fluxes 
if they readjusted target and horn positions using motorize mounts [77J. At NuTeV |216j . 
the monitors demonstrated the alignment of the neutrino beam, necessary for the desired 
precision in siv? 6w The NuMI secondary beam monitors detected misalignment of the 
target which, if uncorrected, could have been catastrophic: the proton beam was initially 
aimed toward the target's edge |219j . which can cause stress on the target fins. The NuMI 
chambers also diagnosed a failure in the target which caused it to fill with cooling water 
|138j . For long-baseline experiments, the muon monitors are essential to demonstrate that 
the neutrino beam points to the remote neutrino detector within the required accuracy, as 
demonstrated at KEK |15j. 

The MiniBooNE beam has a spectrometer system, the "Little Muon Counter" (LMC), 
placed 7° off-axis to its decay pipe (see Figure EZj)- The spectrometer has an upstream 
collimator to define the 7° angle of muons exiting the decay volume, and a small tracking 
system to measure muon momenta exiting the collimator. The experiment, which is designed 
to search for — > oscillations, benefits from direct measurements of intrinsic z/g content in 
the beam, such as originates from /z decays or i^es decays. While the former are constrained 
by measurements of the flux in the MiniBooNE detector, the K^z decays require separate 
knowledge about the i^/vr ratio off the MiniBooNE target. Figure EHl shows the result of a 
Monte Carlo calculation |12(Jj , in the absence of backgrounds from neutrons, conversions, or 
resolution effects from scattering in the collimator, of what would be expected in the LMC 
spectrometer. 
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Figure 49: View of the Fermilab Neutrino Line in 1983, taken from [HZj. The dichromatic 
beam is directed at detectors at 700 m, 1100 m, and 1500 m. 

8 Two-Detector Experiments 

Two-detector experiments were pioneered at Fermilab (see Figure and CERN (see Fig- 
ure EOI) for the purpose of studying neutrino oscillations. Downstream of the neutrino beam, 
a "near" detector measures directly the energy spectrum of neutrinos from the beam. A 
second detector measures the energy spectrum of neutrinos which have propagated for some 
time interval. The distance to the second detector is presumably long compared to the time 
to arrive at the first detector. Deviations between the two energy spectra may be used to 
infer the presence of neutrino oscillations, which manifest themselves as the disappearance 
of the Vy, beam |1961 1331 IT31 1161j . or the appearance of a different neutrino flavor in the z/^ 
beam (eg: Jl]). The direct measurement of the flux in the first detector greatly reduces the 
need to calculate the beam spectrum, improving the experiments' sensitivity. More recently, 
long-baseline experiments have searched for oscillations across distances of order lOO's of km, 
as pioneered at KEK; such has required use of accurate GPS to locate the two separated 
detectors |169j . 

8.1 Calculating the Extrapolated Beam Flux 

Even in the absence of oscillations, the energy spectra in the two detectors are are different, 
and the differences must be calculated so that a "near-to-far" extrapolation of that spectrum 
observed in the near detector can be computed. A comparison of the expected energy spectra 
in the KEK beam at 2;near = 300 m ( "near detector" ) and at zt^^-^ = 250 km ( "far detector" ) 
is shown in Figure |^ The spectra are different, not simply scaleable by the ratio of the 
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Figure 50: Plan view of the CERN PS neutrino beam two-detector experiments, taken 
from [HHl- A bare-target beam taken from the PS is dehvered to the CHARM and CDHS 
experiments. 



two sohd angles subtended by the detectors, -^near/^fai-- What follows is a description of the 
near-to-far extrapolation technique. 

If the beam were a point source, then the prediction could be estimated by 

^far ^ '^FN^ueai 

where the extrapolation factor IZpN = ^nearl'^'jar is just the ratio of solid angles subtended 
by the two detectors. Considering that the beam is an extended source, one could weight 
this extrapolation factor by the pion lifetime along the length of the decay tunnel: 

0.43m7r z 

K-FN — 0.43m^z I^-Jj 

where the integral is over the length L of the decay tunnel and the substitution E^^ ~ Ey/{)A?) 
has been made (c.f. Equation 1^^. The fact that Zp » Z^v reduces the integral in the 
numerator to ~ l/^^, like a point source, while the integral in the denominator reflects the 
more complicated "line source" of neutrinos seen by the near detector. 

Equation 1221 is a simplification, since not all decaying pions produce neutrinos within the 
finite acceptances of the detectors and not all pions are able to decay before interacting along 
the decay pipe walls. Furthermore, the significant acceptance differences between the near 
and far detectors are a function of the energy of the pion (hence neutrino), as is indicated 
schematically in Figure 1^21 Fast pions tend to live longer and decay downstream in the 
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Figure 51: Energy spectra of neutrinos from the KEK beam aX z = 300 m and at 2; = 250 km. 
Taken from [15j. 

decay pipe, closer to the near detector. A variety of pion decay angles 9 will result in a 
neutrino which strikes the transversely large near detector, while only ^ ^ 0° decays result 
in a neutrino arriving at the far detector, affecting the neutrino energy at each detector in 
Equation 1221 Additionally, pions of different momenta enter the decay volume with different 
angular divergences. For an unfocused beam, such as used in the CERN experiments, the 
pion angle leaving the target is 6' ~ 2/7 = 2771.^/ E^. Such differences in angles of entry 
in the decay pipe result (for low-energy pions) in wider-angle decays to reach the neutrino 
detectors, and also greater likelihood of striking the decay volume walls. Even for a horn- 
focused beam, pions of different momenta will have varying divergences depending upon the 
exact tune of the focusing system, and some pions will be unfocused due to the zero-field 
inner-apertures ( "necks" ) of the horns, as shown in Figure |^ 

The calculation of T^pj^ requires a Monte Carlo, as shown in Figure ESI for the NuMI 
beam. The Monte Carlo prediction departs sharply from the idealized curve of Equation I23| 
particularly for very high energy neutrinos (which come from unfocused parent mesons). 




Decay Pipe ND 



Figure 52: Demonstration of the solid angle differences in a two-detector neutrino experi- 
ment. Not to scale. 
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Figure 53: Far-over-near ratio for the NuMI LE beam and the two detector locations at 
z = 1040 m and z = 735.4 km. The points are calculated with a Monte Carlo simulation 
of the beam line, the box sizes correspond to the uncertainties in the calculation due to the 
hadron production model assumed, and the curve is the idealized calculation using only the 
pion lifetime of Equation 1221 



but also at two other distinct values near ^ 4.5 GeV and ^ 8 GeV. Inspection of 
Figure 1^ reveals that these energies are the transition points at which pions no longer pass 
through horn 1 or horn 2, respectively, of the NuMI focusing system. Changes in focusing, 
or the absence of focusing, causes the departure from the idealized curve. Uncertainties 
in the extrapolated far detector flux are typically at the level of (2 — 5)%, well below the 
uncertainty in direct prediction of the flux (which range from (20 — 100)% just from hadron 
production alone), as discussed in Section IH^ Thus, the two-detector experiment serves to 
reduce the uncertainty in the prediction of the "far detector's" energy spectrum. 

Two proposals have been developed to make the far and near detector spectra more 
similar, requiring less sophisticated calculation of T^p^- The T2K collaboration jl31j . noting 
that the structure in T^p^ results from the large angular acceptance of the near detector, 
proposed placing a near detector further from their beam line, at 2; = 2000 m, large compared 
to the 130 m long beam line. The NuMI hadron hose proposal |127j (see Section 16. 2j) was 
developed in part as a means to remove the angular correlations between the near and far 
detectors, thereby the energy difference resulting from Equation |221 The spiraling orbits 
in the hose field randomize the decay angles to the 2 detectors, as well as keep particles 
away from the decay pipe walls so they can continue to follow the pion lifetime curve of 
Equation 1221 [136^. The results from these two ideas are shown in Figures EH and ESI 
respectively. Both proposals are quite effective in removing some of the complexity of the 
near-to-far extrapolation calculation. 
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Figure 54: F/N ratio for the Tokai-to-Kamiokande experiment in Japan. The "far detector" 
is Super Kamiokande at 2; = 250 km. The "near detector" is planned to be at z = 300 m, 
though a proposal has been made to place an additional "near detector" at z = 2000 m. 
Taken from |131j . 
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Figure 55: F/N for the NuMI beam with the addition of the proposed "hadron hose" focusing 
in the decay tunnel. The curve is the idealized calculation using only the pion lifetime of 
Equation |2S1 The hadron hose draws pions away from the decay pipe walls, allowing them 
to decay and thereby following the pion lifetime curve of Equation 1221 Taken from 
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Figure 56: The difference in the energy of the neutrino that a given pion would emit toward 
the far detector to that it would emit toward the near detector in NuMI/MINOS in the low- 
energy (left) and high-energy (right) beams, normalized to the near detector energy. Curves 
are shown with and without the focusing of the Hadron Hose proposal. Taken from |127j . 



That the decay angles to the ND and FD are substantially different is demonstrated in 
Figure In this figure is shown, for each pion focused in the NuMI beam, the ratio of the 
neutrino energy created by that pion if the neutrino strikes the center of the near detector 
at z = 1040 m to the energy of the neutrino created by that same pion if it strikes the 
far detector aX z = 735.4 km. Because of the wide array of decay angles which can strike 
the near detector, particularly for high-energy particles decaying at the downstream end of 
the decay pipe, the energy spectrum in the near detector is systematically lower in the ND 
than the FD. With the inclusion of the Hadron Hose, the ratio of energy spectra at the two 
detectors becomes more symmetric. 

The systematic difference in decay angles motivated another technique jl97j employed for 
analysis of data from the NuMI beam |161j which permits some correction for the near-far 
difference in energy spectra. Noting that a pion decay which gives rise to an energy E^^ 
in the ND actually gives rise to a neutrino energy Ej^ in the FD, a matrix was developed 
which describes the relationship between the energies Ef^^ and Ej^: 

ct^'^^ = S,M,,0r (24) 

where is the flux of neutrinos in the jth energy bin in the FD and 0f ^ is the flux of 
neutrinos in the zth energy bin in the ND. Equation provides a more accurate prediction 
of the flux of neutrinos in the FD given a measurement of the in the ND. Conceptually, 
the matrix is similar to a point-spread function: a particular flux of neutrinos bin of energy 
i in the ND contributes to several bins j in the FD. 
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Neutrino Energy in Near Detector (GeV) 

Figure 57: Graphic representation of the matrix method Mij for predicting the FD flux from 
the ND flux of |197j . The 2-dimensional color plot gives the values of Mij in Equation I24[ 
the scale for which is at the right. The histograms show energy spectra in the ND and FD. 
Figure courtesy M. Messier. 

The matrix elements Mij is shown graphically in Figure EH As noted in the scale, the 
values of Mij are typically of order {zmar / Zfar)"^ = (1/735.4)^. The histograms show energy 
spectra in the ND and FD. As expected, the matrix departs from a diagonal matrix at 
large neutrino energies, which arise from high-momentum pions decaying closer to the ND. 
Summation of the rows of Mij {i.e. summation in j) would give the zth element of the 
far-over- near ratio "^p^- 

Also shown in Figure WJ\ are the energy spectra calculated in the ND and FD. For illus- 
tration, a set of neutrino energies are highlighted in the ND spectrum, and the corresponding 
neutrino energies from the same meson decays are indicated in the FD. A single neutrino en- 
ergy measured in the ND must, because of decay kinematics, contribute to the flux prediction 
for several energy bins in the the FD. Such is the advantage of the matrix method. 
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Figure 58: Density plot showing the pt and Pz of pions which yield neutrinos in the MINOS 
near detector. The box size show the probability of producing such a pion in a p+C collision 
at 120 GeV/c multiplied by the probability for it to go through the beam line and result 
in a neutrino interaction in the MINOS detector. Upper plot: NuMl LE beam; Lower plot: 
NuMl HE beam (see Eigure The symbols indicate previous hadron production mea- 
surements by BartonjlSI; Atherton [221; and NA56/SPY |22]- Forthcoming measurements 
from NA49pT] and MIPP |16Uj will cover much more of the {xf,Pt) plane. 

8.2 Systematic Uncertainties 

This section is not intended to be a full list of all systematic uncertainties for all two- 
detector experiments, nor of all possible techniques which can be used to limit them. Surely 
these will change over time. Rather, the figures shown here are meant to demonstrate where 
systematics appear in the spectrum and what kind of information is helpful in limiting them. 
The NuMl beam will be used as a basis for discussion. 

8.2.1 Hadron Production Uncertainties 

Hadron production data may not exist for a neutrino experiment's target material, thickness, 
beam energy, or relevant portion of {xf,Pt) phase space. Prior to 2006, such was the case 
for NuMl (two-interaction length Carbon target, struck by po=120 GeV/c protons). 

Eigure EHl shows the {xf,Pt) of tt"*" which contribute to the NuMl z/^ energy spectrum 
in either the LE or HE beam configurations. As can be seen, the LE beam configuration 
focuses pions of Pz ~10 GeV/c with px ~ 250 MeV/c. The tail at higher p^ and low pt 
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Figure 59: Energy spectrum of neutrinos interacting on Iron a± z = 735.4 km in the NuMI 
beam calculated in the LE beam configuration (left plot) and the HE beam configuration 
(right plot). Hadron production models from Geant3/Fluka |114j . Malensek |154j . BMPT 
[59] . and MARS-v.14 |157j are compared. 



comes from unfocused pions which pass through the field-free necks of the horns. The HE 
beam focuses pions of ~30 GeV/c. 

In the figure is also shown the data points acquired by high-energy hadron production 
experiments. Only Barton's data is on Carbon, and none are at the correct beam momentum 
or are a thick target measurement. Less constraint is available for HE focusing region, and 
this must be entrusted to hadron production models to extrapolate. No data is available for 
the high energy tail of either of these beams at pt ~ 0. Inspection of Figure EHl shows that 
these regions with little experimental constraint have sizeable model dependence of the fiux 
calculations. The variation amongst various models' predictions is ~20% in the LE beam 
and ~ 30% in the HE beam. 

For a two detector experiment, the uncertainty of relevance is in the far/near ratio, since 
the near detector fiux is measured directly and is used to estimate the far detector fiux. 
Further, the F/N uncertainties, in general smaller than for the direct fiux calculation, are 
only large near the "edges" of the focusing system. If all particles were perfectly focused, 
then hadron production uncertainties would amount to only a rate uncertainty in either 
detector. It is the residual divergence of the beam that results in near-far differences, and 
such arise only at the limits of the focusing or in the high energy tail, where there is no 
focusing. Thus, in Figure ^21 the largest far/near uncertainties (indicated by the vertical 
sizes of the boxes) appear at = 6 GeV, = S GeV, or above E^, > 12 GeV. Beamlines 
with complete focusing for all energies would therefore be anticipated to have reduced model 
dependence of the F/N calculation, such as the Hadron Hose in Figure ESI 
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Figure 60: Distribution of pt and of tt"*" from the NuMI target that contribute to the the 
charged-current event rate at the first detector at z = 1040 m. The box sizes are proportional 
to the probabihty of the pion resulting in a CC interaction in the ND. The 6 plots correspond 
to the 6 beam configurations which have been run: LE10/170kA (top left), LE10/185kA 
(top middle), LE10/200kA (top right), LE100/200kA (bottom left), LE250/200kA (bottom 
middle), LElO/OkA (bottom right) - where each configuration is designated by the number 
of centimeters by which the target is upstream of the horn and by the current in the horns. 
As is evident, each beam configuration samples different region of {xf,Pt)- Taken from |178j . 

For NuMI, it was possible to tune the Monte Carlo hadron production model in [xfiPt) 
to better agree with the near detector energy spectrum by virtue of simultaneously fitting 
data accumulated in several beam energy configurations |161j . lowering the uncertainty in 
the prediction of the far detector spectrum. Similar tuning was used by neutrino experiments 
at BNL [Mi, NOMAD |2H] and NuTeV [122112111, but the NuMI capability to vary target 
position as well as horn current provides additional information. In brief, variation of the 
horn current changes the px kick received by particles through the horn (hence the {pt) of 
the focusing), while the target position dictates the mean xp being focused, as discussed 
in Section m Several focusing uncertainties depend only on p^ oi xp of focused particles, 
as discussed in the next section, so the ability to control both x p and pp helps disentangle 
hadron production and other focusing effects. Figure EOl shows the {pT,xp) sampled by 
several NuMI beam configurations. 
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Figure 61: Uncertainties in tlie near detector spectrum (left plot) and in the far/near ratio 
(right plot) from focusing effects, estimated for the NuMI LE beam. Taken from |178j . 

8.2.2 Focusing Uncertainties 

Several common systematic uncertainties manifest themselves in the detector flux prediction 
or the far/near calculation. Figure 1^ shows some of the larger uncertainties for the NuMI 
LE beam configuration: 

Number of protons on target: two precision toroids, each calibrated by precision 
current sources, were tracked throughout the run against each other and the toroid reading 
the current in the accelerator. A precision of <2% was achieved. 

Proton beam halo scraping on upstream collimating baffles can yield pions from a 
different target location than the nominal target location, and these pions can enter the 
focusing system, contributing high-energy neutrinos. Proton halo was measured by primary 
beam instrumentation to be < 0.2% in magnitude. 

The absolute current in the horns was calibrated by precision toroids around the 
stripline to ±0.5%. 

The current pulse in the horn is ~ 1 msec, at which the skin depth in Aluminum is 
5 = 7 mm, larger than the 3 mm thickness of the horn conductor. The distribution of 
the current in the conductors is therefore not a simple exponential, and the uncertainty on 
this distribution affects strongly those pions which graze and have considerable pathlength 
through the conductor material. 

Misalignments of the horns result in smaller-angle pions, not normally intercepted by 
the horns, receiving some focusing. Such changes the spectrum for those pions at exactly 
the angles near the horn necks. 
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9 Summary 



The present article summarized elements of particle production and focusing to derive con- 
ventional neutrino beams. With 40 years of innovative advances, neutrino beams have gone 
from exploratory devices whose flux and composition was known little better than natural 
neutrino sources to instruments able to predict neutrino fluxes at the experiments at the few 
percent level. Further, the power of neutrino beams has risen from a few detected neutrino 
interactions per week to a few detected interactions per second. 

In the coming years, new facilities will push the technological challenges of neutrino beams 
further. Studies of neutrino oscillations across long baselines of hundreds of kilometers will 
be conducted at the CERN CNGS, Fermilab NuMI, and JPARC-nu facilities. Probing rare 
transitions such as —>■ z/g oscillations, CP violation in the lepton sector, or the appearance 
of the Ur in a t/^ beam will require sophisticated focusing systems capable of coping with 
Megawatts of protons delivered to these facilities over several years. Further, precision 
cross-section experiments being conducted at the Fermilab-MiniBooNE or NuMI beam lines 
or the JPARC-nu beam line will demand accurate demonstration of the neutrino fluxes and 
composition. 

With many exciting fields of physics to be probed using the neutrino, we may look forward 
to continuing advances in beam line systems to realize these goals. 
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A Kinematic Relations 

This appendix reproduces several useful kinematic formulae relating to the energy and decay 
distribution of neutrino daughters from pion and kaon parents. More thorough review of 
relativistic kinematics can be found in |118j . and numerous formulae relevant for neutrino 
beams can be found in |185j . 

Figure 1^21 defines several momentum vectors and directions for the daughters from a vr 
or K decay. Being a two-body decay into a muon and neutrino, the momenta p' of the 
daughters in the center-of-mass frame can be calculated as: 

(25) 

where M is the mass of the it or K parent and is the muon mass. For the daughters in vr 
(K) decay, p' = 29.8 MeV (235.6 MeV). The ir/K being spin zero, the angular distribution 
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Figure 62: View of a parent ir or K meson decay in the laboratory and center-of-mass (CM) 
frames, defining the momenta, energies, and angles of the parent and daughter particles. 



(26) 



of the decay daughters is isotropic in the CM frame, 

dP _ 1 

Transformation of the daughter momenta to the lab is done by the Lorentz boost: 

E = ^iE' + Pp'^) (27) 

P, = l{p', + PE') (28) 

Pt = p't (29) 
where 7 = E/M, = (1 - l/7^)^/^ and E and M are the energy of the parent meson in 
the lab and its mass, respectively. The fact that + = p'^ means that the daughter 
momentum vectors lie on a circle in the CM frame (see Figure . while this relation, upon 
substitution of Equations |2H1 and 1211 yield that the daughter momentum vectors lie on an 
ellipse in the laboratory frame: 

(p, - P^E'^' 



+ ^ = 1 



(30) 



As indicated in Figure ESI the ellipse for the neutrino momentum vector in the lab ap- 
proximately intercepts the origin for very relativistic parents (/5 ~ 1) since E' = p' for the 
neutrino. Because of the large muon energy E'^ (=109 MeV or 258 MeV for ir and K decays, 
respectively), the ellipse for the muon is shifted to the right. As indicated in Figure IHSl the 
neutrino momentum in the lab ranges from to p^"-^ = 2'yp' = (1 — j;p)E (when looking at 
all possible decay angles), which is 0A3E for tt fiu decays and 0.96E for K ^ fiv decays. 
The muon momentum ranges from Q.blE to E imr ^ fiu decays and 0.04ii^ to E in K ^ fiu 
decays. 

Transformation of the daughter angles is found from Equations |2H1 and 1221 noting that 
Pt = P sin 9, p'rp = p' sin 9', p'^ = p' cos 9', and Pz = P cos 9: 

sin 6*' 

7tan6' = , n in,\ (31) 

' coii 9' + {13/ 13') ^ ' 
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CM Frame Laboratory Frame 

Figure 63: Lorentz transformation of a momentum 3- vector p' from the CM frame to the 
laboratory frame. In the CM frame, p' hes on a circle, while in the lab it is required to lie 
on an ellipse (see Equation IHUjl . For (3^1 parents as sketched above, the neutrino ellipse 
approximately is tangent to the pt axis and the muon ellipse is shifted to the right. In 
TT — > /iz/ decays, as sketched above, the muon energy in the lab alwa 



where j3' = p'/E' is the daughter velocity in the CM frame (=1 for the neutrino and =0.28 
or 0.91 for the muon in vr or i^' decays). The maximum decay angle in the lab arises from 
9' = 7r/2, yielding tan6'niax = P' /iP- For the neutrinos from relativistic parents, this reduces 
to 6^^'^ ~ 1/7, for the muons 0^^^ ~ P'/'J- Thus the muons are more forward-boosted. 

The angular distribution of neutrinos in the lab frame is found from Equations [2ni and 
1^ The angular distribution in the lab is found from dP/dVL = [dP / dVL'){dVL' / dVL) = 
{dP/dQ'){d6' /d6) {sin 6' / sin 6). For neutrinos from very relativistic parents with (3^1, 
Equation ED can be inverted to give cos^^' ^ (1 — 7^tan^^)/(l + 7^tan^6'), and 

dP _ 1 472(1 + tan2^)3/2 

dn"^!^ (l + 72tan2^)2 ^ ' 

which reduces to Equation El in the limit that 6 << 1. 

The neutrino energy in the lab is found from Equation EZI For /5 ~ 1, again using 
cos 6' ^ (1 — 7^ tan^ 6)/{l + 7^ tan^ 6) the lab energy is 

1 _ < \ E 

Eu ~ — — 5- — (33) 
1 + 7"^ tan^ 

which reduces to Equation 1^ in the limit that 61, << 1. Equations 1331 and IH^ combine 
to show that the energy distribution of the neutrinos, dP/dEy = {dP / dVLy){dVty / dE^) , is 
constant (which averages over all decay angles O^,). 



69 



References 



[1] CERN Informal Workshop on Neutrino Physics, CERN-63-37, (1963). 

[2] CERN Informal Workshop on Neutrino Physics, CERN-65-32, (1965). 

[3] CERN Informal Workshop on Neutrino Physics, CERN-69-28, (1969). 

[4] l^t International Workshop on Neutrino Beams and Instrumentation, KEK, Tsukuba, 
Japan, 1999, no proceedings published. 

[5] 2^'^ International Workshop on Neutrino Beams and Instrumentation, 6-9 Septem- 
ber, 2000, Fermilab, Batavia, IL, proceedings available online at the following URL: 
[http: //diablo .phys .northwestern, edu/^ngu/biw/program .html , 

[6] 3^"^ International Workshop on Neutrino Beams and Instrumentation, 14-19 March, 
2003, CERN, Geneva, Switzerland, proceedings available online at the following URL: 
[http: //proj-cngs . web . cern. ch/proj-cngs/2002_workshop/announce_l .html 

[7] International Workshop on Neutrino Beams and Instrumentation, 7-11 Novem- 
ber, 2003, KEK, Tsukuba, Japan, proceedings available online at the following URL: 
^ttp : //www-ps . kek. jp/nbi2003/_^ 

[8] S*'' International Workshop on Neutrino Beams and Instrumentation, 7-11 July, 

2005, Fermilab, Batavia, IL, proceedings available online at the following URL: 
^http : / / www . hep . utexas . edu/ nbi2005/ 

[9] 6*^^ International Workshop on Neutrino Beams and Instrumentation, 5-9 September, 

2006, CERN, Geneva, Switzerland, proceedings available online at the following URL: 
thtt p : //pro j -cngs . we b . cern . ch/pro j -cngs/NBI2006/NBI200 6 . htm l 

[10] A.G. Abramov et ai, "Beam Optics and Target Conceptual Design for the NuMI 
Project," Nucl. Instr. Meth. A485, 209 (2002). 

[11] T. Abbot et ai, "Measurement of Particle Production in Proton- Induced Reactions at 
14.6 GeV/c," Phys. Rev. D45, 3906 (1992). 

[12] I. Abt and R. Jongejans, "An absolute calibration of the solid state detectors in the 
narrow band neutrino beam at CERN," Nucl. Instr. Meth. A235, 85 (1985). 

[13] A. Aguilar-Arevalo (MiniBooNE and MINOS Collaborations), "Neutrinos from the 
NuMI beam line in the MiniBooNE detector," talk presented at PANIC '05, Oct. 2005. 

[14] M.H. Ahn et ai, "Search for electron neutrino appearance in a 250 km long baseline 
experiment," Phys. Rev. Lett. 93, 051801 (2004). 



70 



M.H. Ahn et al., "Measurement of Neutrino Oscillation by the K2K Experiment, 
larXiv : hep-ex /0606032[ submitted to Phys. Rev. D (2006). 



L.A. Ahrens et al., "Determination of the Neutrino Fluxes in the Brookhaven Wide 
Band Beams," Phys. Rev. D34, 75 (1986). 

J.V. Allaby et al., "High- Energy Particle Spectra from Proton Interactions at 
19.2 GeV/c," CERN-70-12 (1970). 

J.V. Allaby et al. (CHARM Collab.), "Total Cross- Sections Of Charged Current Neu- 
trino And Anti-Neutrino Interactions On Isoscalar Nuclei," Z. Phys. C38, 403 (1988). 

R.C.Allen et al, Phys. Rev. D47, 11 (1993). 

C. Alt et al, "Inclusive Production of Charged Pions in p + p Collisions at 158 GeV/c 
Beam Momentum," Eur. Phys. J. C45, 343 (2006). 

C. Alt et al., "Inclusive Production of Charged Pions in p+C Collisions at 158 GeV/c 



Beam Momentum," arXiv : hep-ex/0606028 , submitted to Eur. Phys. J. C (2006). 



G. Ambrosini et al., "Measurement of Charged Particle Production from 450 GeV/c 
Protons on Beryllium," Eur. Phys. J. CIO, 605 (1999). 

V.B. Anikeev et al., "Total Cross Section Measurements for v^,V^ Charged Current 
Interactions in 3-30 GeV Energy Range with the IHEP-JINR Neutrino Detector," Zeit. 
f. Phys. C70, 39 (1996). 

D. Antreasyan et al., "Production of vr^ and vr^ at Large Transverse Momentum mp—p 
and p - Colhsions at 200, 300, and 400 GeV," Phys. Rev. Lett. 38, 112 (1977). 

J.G. Asbury et al., "Pion Production at 12° and 15° in Proton-Beryllium Collisions at 
12.5 GeV/c, Phys. Rev. 178, 2086 (1969). 

A. Asner and Ch. Iselin, "A New Focusing System Considerably Increasing the Actual 
PS-Neutrino Beam Flux Intensity," CERN-65-17, May, 1965. 

A. Asner and Ch. Iselin, "Layout of the New CERN Neutrino Beam," CERN-66-24 
(1966). 

P. Astier et al., "Prediction of neutrino fluxes in the NOMAD experiment," Nucl. Instr. 
Meth. A515, 800 (2003). 

H. W. Atherton et al., "Precise Measurements of Particle Production by 400 GeV/c 
Protons on Beryllium Targets," CERN-80-07 (1980). 

B. Aubert et al., "Hadron Production at 0.8 mrad by 300 GeV/c Protons Incident on a 
Thick Aluminum Target," Fermilab-Conf-75/31-Exp (1975). 



71 



[31] W. F. Baker et al, "Particle Production by 10-30 Bev Protons Incident on Al and Be," 
Phys. Rev. Lett. 7, 101 (1961). 

[32] W. F. Baker et al, "Measurement of tt^, p, and p Production by 200 and 300 GeV/c 
Protons," Phys. Lett. B51, 303 (1974). 

[33] W. F. Baker et al, "Measurement of n'^, K^, p, and p Production by 400 GeV/c 
Protons," Fermilab-Conf-78/79-Exp (1978). 

[34] D.G. Baratov et al., "Loading of Parabolic Lenses by a High Electric Current," Zh. 
Tekh. Fiz. 45, 1954 (1974) [trans. Sov. Phys. Tech. Phys. 20(9), 1229 (1976)]. 

[35] D.G. Baratov et al., "Neutrino Beam of the Institute of High Energy Physics I: General 
Description of and Basic Characteristics," Zh. Tekh. Fiz. 47, 991 (1977) [trans. Sov. 
Phys. Tech. Phys. 22(5), 591 (1977)]. 

[36] D.G. Baratov et al, "Neutrino Beam of the Institute of High Energy Physics II: Char- 
acteristics of the 70-GeV Proton Beam Extracted for Neutrino Experiments," Zh. Tekh. 
Fiz. 47, 997 (1977) [trans. Sov. Phys. Tech. Phys. 22(5), 595 (1977)]. 

[37] D.G. Baratov et al, "Neutrino Beam of the Institute of High Energy Physics HI: System 
for Shaping the Proton Beam for the External Target," Zh. Tekh. Fiz. 47, 1007 (1977) 
[trans. Sov. Phys. Tech. Phys. 22(5), 601 (1977)]. 

[38] D.G. Baratov et al, "Neutrino Beam of the Institute of High Energy Physics IV: Focus- 
ing Device of the Neutrino Channel," Zh. Tekh. Fiz. 47, 1014 (1977) [trans. Sov. Phys. 
Tech. Phys. 22(5), 605 (1977)]. 

[39] D.G. Baratov et al, "Neutrino Beam of the Institute of High Energy Physics V: Pulsed 
Power Supply for the Neutrino Focusing System," Zh. Tekh. Fiz. 48, 91 (1978) [trans. 
Sov. Phys. Tech. Phys. 23(1), 53 (1978)]. 

[40] D.G. Baratov et al., "Neutrino Beam of the Institute of High Energy Physics VI: Systems 
for Optimization and Control in Neutrino Experiments," Zh. Tekh. Fiz. 48, 99 (1978) 
[trans. Sov. Phys. Tech. Phys. 23(1), 58 (1978)]. 

[41] D.G. Baratov et al., "Neutrino Beam of the Institute of High Energy Physics VII: 
Particle Loss and Radiation Environment Along the Transport System for the Extracted 
Proton Beam in the Neutrino Channel," Zh. Tekh. Fiz. 48, 109 (1978) [trans. Sov. Phys. 
Tech. Phys. 23(1), 64 (1978)]. 

[42] B.C. Barish et al. (CITF CoUab.), "Exploratory Study of High- Energy Neutrino Inter- 
actions," Phys. Rev. Lett. 31, 565 (1973). 

[43] B.C. Barish et al. (CITF CoUab.), "Neutral Currents in High-Energy Neutrino Colli- 
sions: An Experimental Search," Phys. Rev. Lett. 34, 538 (1975). 



72 



[44] B.C. Barish et al, "Charged Current Neutrino And Anti-Neutrino Cross-Section Results 
From The CITFR Experiment" Fermilab-Conf-78-046-E, Presented at 3rd Int. Conf. on 
New Results in High Energy Physics, Nashville, Tenn., Mar 6-8, 1978. 

[45] B.C. Barish, "Experimental Aspects Of High-Energy Neutrino Physics," Physics Re- 
ports 39, 279 (1978). 

[46] B.C. Barish et al, "Recent Results on Total Neutrino and Antineutrino Cross Sec- 
tions by the CFRR Collaboration," Fermilab-81 /80-EXP, presented at SLAC Summer 
Institute on Particle Physics, 1981. 

[47] S.J. Barish et al, "Study of Neutrino Interactions in Hydrogen and Deuterium: De- 
scription of the Experiment and Study of the Reaction i/ + d ^ ii~ + p + Ps" Phys. 
Rev. D16 3103 (1977). 

[48] D.S. Barton et al, "Experimental Study of the A Dependence of Inclusive Hadron Frag- 
mentation," Phys. Rev. D27, 2580 (1983). 

[49] A.A. Batalov et al, "A focusing Device for Neutrino Sources at Meson Factories," Nucl. 
Instr. Meth. A251, 231 (1986). 

[50] D. Beavis et al "P889: Long Baseline Neutrino Oscillation Experiment at the AGS," 
Report No. BNL-52459, April, 1995. 

[51] A. Benvenuti et ai (HPWF Collab.), "Early Observation of Neutrino and Antineutrino 
Events at High Energies," Phys. Rev. Lett. 30, 1084 (1973). 

[52] A. Benvenuti et al (HPWF Collab.), "Observation of Muonless Neutrino- Induced In- 
elastic Interactions," Phys. Rev. Lett. 32, 800 (1974). 

[53] A. Benvenuti et al (HPWF Collab.), "Observation of New Particle Production by High- 
Energy Neutrinos and Antineutrinos," Phys. Rev. Lett. 34, 419 (1975). 

[54] J. P. Berge et al (CDHS Collab.), "Total Neutrino and Antineutrino Charged Cur- 
rent Cross Section Measurements in 100, 160, and 200 GeV Narrow Band Beams," Z. 
Phys.C35, 443 (1987). 

[55] F. Bergsma et al, "A Search for Neutrino Oscillations," Zeit. f. Phys. C40 171 (1988). 

[56] R. Bernstein et al, "Sign-Selected Quadrupole Train," Fermilab-TM-1884 (1994). 

[57] R. Blair et al, " Monitoring And Calibration System For Neutrino Flux Measurement 
In A High-Energy Dichromatic Beam," Nucl. Instr. Meth. A226, 281 (1984). 

[58] D. Bloess et al, "Determination of the Neutrino Spectrum in the CERN 1967 Neutrino 
Experiment," Nucl. Instr. Meth. 91, 605 (1971). 



73 



[59] M. Bonesini, A. Marchionni, F. Pietropaolo, and T. Tabarelli de Fatis, "On Particle 
Production for High Energy Neutrino Beams," Eur. Phys. J. C20, 13-27 (2001). 

[60] A.E. Brenner et al, "Experimental Study of Single-Particle Inclusive Hadron Scattering 
and Associated Multiplicities," Phys. Rev. D26, 1497 (1982). 

[61] S. Brice et ai, " A Proposal for a near detector experiment on the booster neu- 
trino beamline: FINeSSE: Fermilab intense neutrino scattering scintillator experiment," 
Fermilab-Proposal-937 (2004). 

[62] L. Bruno, "The CNGS Target," in ||. 

[63] G.l. Budker, International Conference on Accelerators, Dubna, August 1963, pp. 282- 
283. 

[64] B.F. Bayanov, G.I. Budker, et ai, "Wide- Angle Optics with Strong Magnetic Fields for 
Efficient Generation of Secondary-Particle Beams," Proc. 10 Int. Conf. High Energy 
Accel., Protvino, Vol.2, pp. 109-109 [trans, available as SLAG preprint SLAG-TRANS- 
0213 (1984)]. 

[65] Bugorsky, A. P. et al, "Muon Flux Measuring System for Neutrino Experiments at the 
IHEP Accelerator," Nucl. Inst. Metk, 146, 367-371 (1977). 

[66] R. Burns et ai, "Search for Intermediate Bosons in High- Energy Neutrino Interactions", 
Phys. Rev. Lett. 15, 42 (1965). 

[67] R. Burns et el., "Determination of Neutrino Flux," in pp. 97-110. 

[68] J. Gamas et ai, "Screens Versus SEM Grids for Single Pass Measurements in SPS, 
LEP, and LHG," GERN SL/95-62 (BI), presented at the Second European Workshop on 
Beam Diagnostics and Instrumentation for Particle Accelerators (DIPAG'95), Lubeck- 
Travemunde, Germany, 28-31 May 1995. 

[69] U. Gamerini et ai, "An Experiment to search for neutrino oscillations using a electron 
- neutrino enriched beam," Fermilab-PROPOSAL-0664, 1980. 

[70] U. Gamerini et ai, "An Experiment to search for muon-neutrino/electron-neutrino 
tau- neutrino neutrino oscillations using an enriched electron- neutrino/electron- 
antineutrino beam," Fermilab-PROPOSAL-0665, 1980. 

[71] J. Gampbell et al, "Study of the Reaction up — * /x^vr+p," Phys. Rev. Lett. 30, 335 
(1973). 

[72] D.G. Garey, L.J. Edwards, Y.W. Kang, F.A. Nezrick, R.J. Stefanski, and T.E. Toohig, 
"500 GeV Neutrino Beam at NAL," Fermilab-TM-472 (1970). 



74 



[73] D.C. Carey, R.J. Stefanski, and L.C. Teng, "Wide Band Neutrino Beams with 
Quadrupole Focusing," IEEE Trans. Nucl. Sci. 18, 755 (1971), also Fermilab-Conf- 
71/045. 

[74] D.C. Carey, "Neutrino Optics," Nucl. Instr. Meth. A363, 162 (1995). 

[75] A. Carroll et al., "Large Acceptance Magnetic Focusing Horns for Production of a High 
Intensity Narrow Band Neutrino Beam at the AGS," IEEE Trans. Nucl. Sci. 32, 3054 
(1985). 

[76] A. Carroll, et ai, "Overview of Recent Focusing Horns for the BNL Neutrino Program," 
Proc. IEEE Part. Accel. Conf., Washington, D.C, p. 1731 (1987). 

[77] Casagrande, L. et ai, "The alignment of the CERN West Area neutrino facility," CERN 
Yellow Report 96-06 (1996). 

[78] M.G. Catanesi et ai, "Measurement of the production cross-section of positive pions in 
pAl collisions at 12.9 GeV/c," Nucl. Phys. B732, 1 (2006). 

[79] G. Cavallari, E.H.M. Heijne, P. Jarron, P. Lazeyras, M. Musso, "Solid state detectors 
used for the CERN Neutrino Flux Monitoring (NFM)" IEEE Trans. Nucl. Sci. 25, 600 

(1978) 

[80] C.Y. Chi et ai, "Segmented Ionization Chambers for High Intensity Beam Measure- 
ments," Nucl. Instr. Meth. A281, 448 (1985). 

[81] Y. Cho et ai, "Pion Production in Proton-Beryllium Collisions at 12.5 GeV/c, Phys. 
Rev. D4, 1967 (1971). 

[82] G. Cocconi, L. Koester, and D.H. Perkins, LBL Report No. UCRL 10022 (1961). 

[83] J. Conrad, M. Shaevitz, T. Bolton, "Precision Measurements with High-Energy Neu- 
trino Beams," Rev. Mod. Phys. 70, 1341 (1998). 

[84] G. Danby et ai, "Observation of High-Energy Neutrino Reactions and the Existence of 
Two Kinds of Neutrinos" , Phys. Rev. Lett. 9, 36 (1962). 

[85] G.T. Danby et al, Bull. Am. Phys. Soc. 10, 34 (1965). 

[86] LA. Danilchenko, R.A. Rzaev, A.V. Samoylov, and V.I. Voronov, "Focusing System for 
Neutrino Parents," Nucl. Instr. Meth. 105, 147 (1972). Also Internal Serpukhov Report 
IHEP-70-93 (1970). 

[87] L.L. Danilov et al, "Parabohc Lens for Axially- Symmetric Focusing of Relativistic 
Particles," Zh. Tekh. Fiz. 37, 914 (1967) [trans. Sov. Phys. Tech. Phys. 12(5), 656 
(1967)]. 



75 



D. Dekkers et al, CERN Internal Report NP/Int./64-5 (May, 1964), later published as 
D. Dekkers et al, Phys. Rev. 137, 962 (1965). 

G. De Lellis, P. Migliozzi and P. Santorelli, "Charm Physics with Neutrinos," Physics 
Reports 399, 227 (2004). 

M. Derrick, "Plans for Neutrino Experiments with the 12 ft. Bubble Chamber at the 
ZGS," in 0, pp. 3-12. 

A. Diddens et ai, "N4 Experiment at CERN," unpublished. 

M. Diwan, S. Kahn, R.B. Palmer, "A Solenoidal Capture System for Neutrino Produc- 
tion," Proa. IEEE Part. Accel. Conf., New York, NY (1999). 

V. Dohm et ai, "Linear Focusing of Neutrino Parents," Nucl. Instr. Meth. 124, 477 
(1975). 

Need a reference for DPMJET-III. 

J.C. Dusseux, J.B.M. Pattison, and G. Ziebarth, "The CERN Magnetic Horn and its 
Remote Handling System," CERN-72-11 (1972). 

D.A. Edwards and F.J. SciuUi, "A Second-Generation Narrow Band Neutrino Beam," 
Fermilab-TM-660 (1976). 

D.A. Edwards and F.J. SciuUi, "350 GeV/c Dichromatic Neutrino Target Train," 
Fermilab-TM-661 (1976). 

R.T. Edwards et al, "Forward Inclusive Production of Ks, A°, A°, and n in the collision 
of 200 GeV/c tt", fsT", p, and p on Be," Phys. Rev. D18, 76 (1978). 

T. Eichten et al, "Particle Production in Proton Interactions in Nuclei at 24 GeV/c" 
Nucl. Phys. B44, 333 (1972). 

[100] T. Eichten et ai, "High Energy Electronic Neutrino (z/g) and Antineutrino (z/g) Inter- 
actions," Phys. Lett. B46, 281 (1973). 

[101] K. Elsener, "The CNGS Beam," in 0. 

[102] G. Ferioli and R. Jung, "Evolution of Secondary Emission Efficiencies of Various Ma- 
terials Measured in the CERN SPS Secondary Beam Lines," CERN-SL-97-71(BI), pub- 
lished in Proceedings of Beam Diagnostics and Instrumentation for Particle Accelerators 
(DIPAC), Frascati (Rome), Italy, Oct. 12-14, 1997. 

[103] J. Simon, J. Crimson, J. Lindberg, D. Theriot, "Design of the Fermilab Remote Target 
Maintenance System," Fermilab-TM-0618 (1975); "Fermilab Target Areas and Target 
Train Systems," Fermilab-TM-0616 (1975); "Design of the Fermilab Neutrino Horn 
Train System Referenced to Radioactive Maintenance," Fermilab-TM-0617 (1975). 

76 



[104] R.P. Feynman, "Very High-Energy Collisions of Hadrons," Phys. Rev. Lett. 23, 1415 
(1969). 

[105] H.E. Fisk and F. Sciulli, "Cliarged-Current Neutrino Interactions," Annual Review of 
Nuclear and Particle Science. 32,499 (1982). 

[106] G. Collazuol, A. Ferrari, A. Guglielmi, P.R. Sala, Nucl. Instr. Meth. 449, 609 (2000). A. 
Ferrari and P.R. Sala, ATLAS Internal Note ATL-PHYS-97-113, Proc. of the Workshop 
on Nuclear Reaction Data and Nuclear Reactors Physics, Design and Safety, ICTP, 
Trieste, Italy 1996. Publ. by World Scientific, A.Gandini, G.Reffo, eds.; P. Sala, talk 
given at the 0. 

[107] E.B. Forsythe, L.M. Lederman, and J. Sunderland, "The Brookhaven- Columbia 
Plasma Lens," in 0, pp. 131-140. 

[108] P. Fritze et al, Phys. Lett. 96B, 427 (1980). 

[109] J.M. Gaillard, "The Brookhaven Neutrino Experiment," in 0, pp. 33-42. 

[110] V. Garkusha, private communication. See also V. Garkusha et al, IHEP-85-82 (1985). 

[Ill] C. Zeitnitz and T.A. Gabriel, Nucl. Instr. Meth. A349, 106 (1994). 

[112] S. Agostinelh et al. (GEANT4 Collaboration), Nucl. Instr. Meth. A506, 250 (2003). 

[113] H. Fesefeldt, PITHA-85-02. 

[114] GEANT Detector Description and Simulation Tool. CERN Program Library, W5013 
(1994). 

[115] M. Giesch et al, "Status of magnetic horn and neutrino beam," Nucl. Instr. Meth., 20, 
58 (1963). 

[116] J. Crimson and S. Mori, "New Single Horn System," Fermilab-TM-824 (1978). 

[117] H. Grote, R. Hagedorn, J. Ranft, "Particle Spectra," CERN (1970). 

[118] R. Hagedorn, Relativistic Kinematics, (New York: John Wiley & Sons), 1963. 

[119] D.A. Harris et al., "Precision calibration of the NuTeV calorimeter," Nucl. lustrum. 
Meth. A447, 377 (2000). 

[120] T. Hart, "The MiniBooNE Little Muon Counter," in |H]. 

[121] E.H.M. Heijne, "Influence of muon induced secondary radiation on the muon flux 
measurement in the CERN neutrino beams," CERN- 79-4 (1979). 



77 



[122; 
[123; 

[124 

[125' 
[126 

[127; 
[128; 



[129 

[130 
[131 

[132 
[133 

[134 

[135 

[136 
[137 



E.H.M. Heijne, P. Jarron, P. Lazeyras, W.R. Nelson, G.R. Stevenson, "A tiny telescope 
of Si-detectors for high energy muon flux measurement with electron rejection," IEEE 
Trans. Nucl. Sci. 27, 272 (1980). 

E.H.M. Heijne, "Muon flux measurement with silicon detectors in the CERN neutrino 
beams," CERN-83-06 (1983). 

Hill, J. E. et al, "Beam monitoring and modeling for the K2K long-baseline neutrino 
oscillation experiment," Int. J. Mod. Phys., A16S1B, 758-760 (2001). 

M. Holder, "Experimental Determination of the Neutrino Flux," in |2I, pp. 115-120. 

S. Humphries, Principles of Charged Particle Acceleration (New York: John Wiley &; 
Sons), 1986. 

J. Hylen et al, "The Hadron Hose: Continuous Toroidal Focusing for Conventional 
Neutrino Beams," Nucl. Instr. Meth. A498, 29 (2003). 

Y. Itow et al, "The JHF-Kamioka Neutrino Project," KEK-REPORT-2001-4, ICRR- 
REPORT-477-2001-7, TRI-PP-01-05, Prepared for 3rd Workshop on Neutrino Os- 
cillations and Their Origin (NOON 2001), Kashiwa, Japan, 5-8 Dec 2001. e-Print 
arXiv:hep-ex/0106019. 



D. E. Jaffe et al., "Comparison of inclusive particle production in 14.6-GeV/c proton- 
nucleus collisions with simulation," Nucl. lustrum. Meth. B246, 309 (2006) 

M. Jonker et al, Phys. Lett. 96B, 435 (1980). 

E. Kearns et al, "A Proposal for a Detector 2 km Away From the T2K Neutrino 
Source," presented to the Neutrino Scientific Assessment Group subpanel of U.S. Dept. 
of Energy High Energy Physics Advisory Panel, May 30, 2005. 

T. Kirk, "BNL E910," in 0. 

T. Kitagaki et al, "Neutrino Flux and Total Charged- Current Cross Sections in High- 
Energy Neutrino-Deuterium Interactions," Phys. Rev. Lett. 49, 98 (1982). 

K. Kodama et al, "Observation of tau neutrino interactions," Phys. Lett. B504 218 
(2001). 

D.D. Koetke et al, "A Direct-Current Pion Focusing Magnet for Low-Energy In-Flight 
Muon Neutrino Beams," Nucl. Instr. Meth. A378, 27 (1996). 

S. Kopp, "The NuMI Hadron Hose," presented at jS]. 

S. Kopp et al, "Segmented Foil SEM Grids at Fermilab," Fermilab-Conf-05-092, sub- 
mitted to 2005 IEEE Nuclear Science Symposium, Puerto Rico. 



78 



[138 

[139 

[140 
[141 

[142; 

[143 

[144 

[145 

[146 

[147 

[148 

[149 

[150 

[151 

[152; 

[153 



S. Kopp et al, "Secondary Beam Monitors for the NuMI Facility at FNAL," Fermilab- 
Pub-06-007-AD (2006), accepted in Nucl. Instr. Meth. A. 

M. Kostin, S. Kopp, M. Messier, D. Harris, J. Hylen, A. Para, Fermilab-TM-2353-AD 
(2002), "Proposal for Continuosly- Variable Beam Energy". 

I. Kourbanis, "MiniBooNE Horn and Target System," in [Bj. 

R.L. Kustom et al., "Quasielastic Neutrino Scattering," Phys. Rev. Lett. 22, 1014 
(1969). 

"Proc. Workshop on Prog. Options in Interm. Energy Physics," J.C Aired and B. 
Talley, eds., Los Alamos Conf. Proc. LA-8335-C (Los Alamos, 1979). 

"Proc. of the Workshop on Nucl. and Part. Phys. up to 31 GeV," J.D. Bowman, L.S. 
Kisslinger, and R.R. Silbar, Los Alamos Conf. Proc. LA-8775-C (Los Alamos, 1981). 

"Proc. of the Los Alamos Neutrino Workshop," F. Boehm and G.J. Stephensen, Jr., 
eds., Los Alamos Conf. Proc. LA-8775-C (Los Alamos, 1981). 

"A Proposal to the Department of Energy for a High Intensity Los Alamos Neutrino 
Source from the Los Alamos National Laboratory," (Los Alamos, 1982) (unpublished). 

L. Lederman, "Neutrino Experiments at Brookhaven National Laboratory," in |2|, pp. 
1-10. 

P. Limon et al, "A Sign-Selected Di-Chromatic Neutrino Beam," Nucl. Instr. Meth. 
116, 317 (1974). 

R.L. Burman & W.C. Louis, "Neutrino Physics at Meson Factories and Spallation 
Neutron Sources," J. Phys. G29, 2499 (2499). 

R.A. Lundy et al, "71"*= and Production Cross Sections for 12.5 BeV Protons on 
Be," Phys. Rev. Lett. 14, 504 (1965). 

J. Link, "Pion Production on Be at 6, 12, 18 GeV at E910," talk presented at NuFact04, 
Osaka, Japan (2004). 

K.T. McDonald, "A Neutrino Horn Based on a Solenoid Lens," 
larXiv: physics/0312022 (2003). 

K. McFarland et al, "Measurement of sin^ ^vi/ from Neutrino Nucleon Scattering at 
NuTeV," 33rd Rencontres de Moriond: Electroweak Interactions and Unified Theories, 
Les Arcs, France, 14-21 Mar 1998,| arXiv:he p-ex/9806013l 

R.D. McKeown and P. Vogel, "Neutrino Masses and Oscillations: Triumphs and Chal- 
lenges," Physics Reports, 394, 315 (2004). 



79 



[154] A.J. Malensek, "Empirical Formula for Thick Target Particle Production," Fermilab- 
FN-341 (1981). 



[155; 
[156 

[157; 



[158 
[159 
[160 
[161 

[162 
[163 

[164 

[165 

[166 
[167; 

[168 

[169 

[170 



A. Mann et al, "Study of the Reaction un /i-p," Phys. Rev. Lett. 31, 884 (1973). 

G. J. Marmer et al., "Differential Production Cross Sections of Low-Momentum Parti- 
cles from 12.3 BeV/c Protons on Beryllium and Copper," Phys. Rev. 179, 1294 (1969). 

N.V. Mokhov, "The MARS Monte Carlo", Fermilab FN-628 (1995); O.K. Krivosheev 
et al, Proc. of the Third and Fourth Workshops on Simulating Accelerator Radiation 
Environments (SARE3 and SARE4), Fermilab-Conf-98/043(1998) and Fermilab- Conf- 
98/379(1998). 

T. Maruyama, "Measurement of Neutrino Spectrum by Pion Monitor at K2K," in [S]. 
M. Messier, "Monte Carlo Simulations of the NuMI Beam," in jS]. 
M. Messier, "The MIPP/E907 Experiment at Fermilab," in |H]. 

D.G. Michael et al, "Observation of Disappearance of Muon Neutrinos with an Accel- 
erator Neutrino Beam," Fermilab-Pub-06-243, submitted to Phys. Rev. Lett. (2006). 

R.H. Milburn, "Theory of the Hadronic Hose," Fermilab-NuMI-B-271 (2000). 

P. Musset and J. -P. Vialle, "Neutrino Physics with Gargamelle," Physics Reports 39, 
1-130 (1978). 

S. Mori, "Wide-Band Single Horn System," Fermilab-TM-663 (1976); "Wide-Band 
Single Horn System (II)," Fermilab-TM-720 (1977). 

S. Mori, S. Pruess, and R. Stefanski, "Electron Neutrino Beam," Fermilab-TM-725 
(1977). 

S. Mori "Improved Electron Neutrino Beam," Fermilab-TM-769 (1979). 

F.A. Nezrick, "A monoenergetic neutrino beam using current-sheet focusing elements," 
IEEE Trans. Nucl. Sci. 18, 759 (1971). 

F.A. Nezrick, "Fermilab Neutrino Horn Focusing System," IEEE Trans. Nucl. Sci. 22, 
1479 (1975). 

H. Noumi et al., " Precision positioning of Super-Kamiokande with GPS for a long 
baseline neutrino oscillation experiment," Nucl. lustrum. Meth. A398, 399 (1997). 

D. Ayres et al., "Letter of Intent to build an Off- axis Detector to study numu to nue 
oscillations with the NuMI Neutrino Beam," ,arXiv:hep-ex/0210005 . 



80 



[171 

[172 

[173; 
[174 

[175; 

[176 

[177 

[178 

[179 
[180 

[i8i; 

[182 
[183 
[184 

[185 

[186 
[187 
[188 



R.A. Ong, L. Sulak, R. Watts, "Design of an Enhanced 1 GeV Electron Neutrino 
Beam," Nucl. Instr. Meth. A236, 256 (1985). 

A. Orkin-Lecourtois et ai, "Fluxes of Neutrinos and Antineutrinos Produced at 
CERN," in [2|, pp. 111-114. 

R.B Palmer, "Magnetic Fingers," in 0, pp. 141-146, Jan. 1965. 

W.K.H. Panofsky and W.R. Baker, "A Focusing Device for the External 350 MeV 
Proton Beam of the 184-inch Cyclotron at Berkeley," Rev. Sci. Instr. 20, 445 (1950). 

W.K.H. Panofsky, "Remarks on the 'Horn of Plenty' for Neutrino Experiments," 
SLAC-TN-62-11 (1962). 

A. Para and M. Szleper, "Neutrino Oscillations Experiments Using Off-axis NuMI 
Beam," FERMILAB-PUB-01-324 (2001) arXiv:hep-ex/0 1100321 

J.M. Pattison, "The CERN Neutrino Beam," in pp. 13-32. 

Z. Pavlovic, "Studies of the Neutrino Flux for the NuMI Beam at FNAL," PhD Thesis, 
University of Texas at Austin. 

G. Plass and B. Kuiper, "Fast Beam Extraction and the External Proton Beam," in 
[IJ, pp. 83-96. 

G. Plass, K.M. Vahlbruch, and H.W. Wachsmuth, "The Determination of the Spectra 
of High- Energy Neutrino Parents," in j2], pp. 121-128. 

G. Plass and H.W. Wachsmuth, "Muon Flux Measurement in the Neutrino Shielding," 
in 0, pp. 147-152. 

G. Plass, K.M. Vahlbruch, and H.W. Wachsmuth, "Test of a Spectrometer for Intense 
Nanosecond Bursts of High- Energy Particles," Nucl. Instr. Meth. 46, 223 (1967). 

B. Pontecorvo, J. Exp. Theor. Phys. (U.S.S.R.) 37, 1751 (1959) [tranlation: Sov. 
Phys.-JETP 10, 1236 (I960)]. 

B. Pontecorvo, " 'Direct' Neutrino Production and Charmed Particles," Zh. Eksp. 
Teor. Fiz. 69, 452 (1975) [trans. Sov. Phys.-JETP 42(2), 229 (1975)]. 

C. Ramm, "Neutrino Spectra from the Two-Body Decay of Relativistic Parents," in 
[H, pp. 111-131. 

S. Rangod, "Horn Design for the Neurtino Factory," in [H]. 

E. Regenstreif, "Contribution to the Theory of the Beam Guide," CERN-64-41 (1964). 
A. Roberts, "Simple High-Momentum Neutrino Beams," Fermilab-FN-124 (1968). 



81 



[189] J.R. Sanford and C.L. Wang, AGS Internal Reports, BNL-11479 & BNL-11299, 
Brookhaven National Laboratory, 1967 (unpublished). 



[190 
[191 
[192 

[193 

[194 
[195 

[196 

[197; 

[198 
[199 

[200 

[201 
[202 
[203 

[204 
[205 
[206 
[207; 



D. Schmitz, "HARP Data on the MiniBooNE Target," in (Hj. 
M. Schwartz, Phys. Rev. Lett. 4, 306 (1960). 

W.P. Sims et ai, "Assembly Techniques Used in Construction of Neutrino Horns at 
BNL," Proc. IEEE Part. Accel. Conf., p. 1734, Washington, DC, 1987. 

P. Skubic et al, "Neutral-Strange-Particle Production by 300 GeV/c Protons," Phys. 
Rev. D18, 3115 (1978). 

R. Stefanski and H. White, "Neutrino Flux Distributions," Fermilab-FN-292 (1976). 

J. Steinberger, "Experiments with High-Energy Neutrino Beams," Nobel Lecture, De- 
cember, 1988. 

I. Stockdale et ai, "Limits on Neutrino Oscillations in the Mass Range 55 < Am^ < 
800 eV^," Phys. Rev. Lett. 52, 1384 (1984). 

M. Szleper and A. Para, "Neutrino Spectrum at the Far Detector Systematic Errors," 
Fermilab-Pub-0 1-309 (2001). 

G. Tassotto, "Fermilab Primary Beam Instrumentation," in 0. 

S. van der Meer, "Comparison between quadrupole lenses of classical design and 'Panof- 
sky lenses'," CERN-PS-INT-MM-59-8 (1959). 

S. van der Meer, "A directive device for charged particles and its use in an enhanced 
neutrino beam," CERN-61-07 (1961). 

S. van der Meer, "The Beam Guide," CERN-62-16 (1962). 

S. van der Meer, "Neutrino Flux Calculations," in [T], pp. 97-110. 

S. van der Meer, "Improved Collection of Secondaries from a Long, Narrow Target by 
a Horn Doublet," CERN/PS/AA/80-12 (1980). 

W.A. Venus, "Scaling the CERN Neutrino Beam to Higher Energies," in P], pp. 59-66. 

H. Vogel et ai, Proc. Int. Symp. on Magnet Tech., SLAC (1965), pp. 650-6. 

G. von Dardel, CERN-NP/Internal Report 62-17. 

H. Wachsmuth, "Neutrino Flux Measurement," in p. 41. 



82 



[208] H. Wachsmuth, "400 GeV Proton-Induced Prompt Neutrinos at and 15 Milliradian," 
presented at Topical Conference on Neutrino Physics, Oxford, CERN/EP/PHYS 78-29 
(1978). 

[209] H. Wachsmuth, "Past and Future Oscillation Experiments in CERN Neutrino Beams," 
Proc. Neutrino 1982 Conference, Belatonfured, Germany, CERN/EP-82-102 (1982). 

[210] J.C. Walker et al., "Mechanical Support and Transport System Used for the Neutrino 
Horn System at Brookhaven National Laboratory," Proc. IEEE Part. Accel. Conf., 
Washington, D.C., p. 1737 (1987). 

[211] C.L. Wang, "Pion, Kaon, and Antiproton Production Between 10 and 70 BeV," Phys. 
Rev. Lett. 25, 1068 (1970); erratum Phys. Rev. Lett. 25, 1536 (1970). 

[212] C.L. Wang, "Empirical Formula for Pion Production in Proton-Proton Collisions up 
to 1500 GeV," Phys. Rev. D7, 2609 (1973). 

[213] C.L. Wang, "Pion Production in High-Energy Collisions," Phys. Rev. DIO, 3876 
(1974). 

[214] H. White, "Double Horn Neutrino Production with a Test Target," Fermilab-TM-662 
(1976). 

[215] Y. Yamanoi, "Large Horn Magnets at the KEK Neutrino Beam Line," KEK-Preprint- 
97-225, presented at 15*^ Int. Conf. on Magnet Tech., Beijing, China, Oct. 1997. 

[216] J. Yu et al, "NuTeV SSQT Performance," Fermilab-TM-2040 (1998). 

[217] G.P. Zeller et al, "A Precise Determination of Electroweak Parameters in Neutrino 
Nucleon Scattering," Phys. Rev. Lett. 88, 091802 (2002). 

[218] R. Zwaska, "Accelerator Systems and Instrumentation for the NuMI Beam at FNAL," 
PhD Thesis, University of Texas at Austin, December, 2005. 

[219] R. Zwaska, et al, "Beam-Based Alignment of the Target Station Components of the 
NuMI Facihty at FNAL," Fermilab-Pub-06-171-AD (2006), submitted to Nucl. Instr. 
Meth. 



83 



Problem Set 



1. Assume a "half-horn" with an inner-conductor which is partially tapered and at z = / is 
perpendicular to the z axis, as in FigurelMl We will calculate the trajectory of a particle 
of momentum p through such a horn but not resort to the "thin lens" approximation 
that the pathlength s through the horn is small compared to the distance Zin from the 
particle source (target) to the horn. 

(a) Show that, if the particle enters the horn at a coordinate {^'y^ir-y^^ its trajectory 
through the horn is given by 

r = rj^exp[ylo(cos^^ — cos6'jj^)] 
^ ~ ^in = ^oHn / (i6' cos 6'exp[y4o (cos 6^ — cos ft^))] 

J6»=0 

where rj^/^j^ = tan^j^^ and = = p(GeV/c)/(6 x 10-5/(kA)). The above 
are the equations integrated numerically for several values of ^j^^ in Figure El 

(b) If the particle angles entering this horn are small, namely << 1 and ^o^j^ << 
1, show that the entrance points of the particles into the horn should follow: 

(-in - \ff + = 

which is an ellipse with major and minor axes given by a = |/, b = 2a/(3Ao)^ . 

(c) Show that if instead one requires a tighter "small-angle" approximation ^o^f^^ << 
1, one obtains the condition of a parabola. 



What is the coefficient A? The above form is for a parabolic lens. 
r \ 




Figure 64: Schematic of a particle trajectory in a lens not assumed to be "thin." 1. -upstream, 
tapered end of inner conductor, 2. -downstream end of inner conductor, which points per- 
pendicular to the z axis, 3.-source of particles (target) at {z,r) = (0, 0). 
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Figure 65: End-view of a quadrupole magnet of aperture a. 

2. This problem concerns the wide-band focusing properties of quadrupole magnets. 

(a) Show that the focal length of a quadrupole of length L, aperture a, and maximum 
field Bq (in Tesla) at the pole tip, acting on a momentum p (in GeV/c) is 

/ 



O.SBqL 

The magnetic potential of a quadrupole is conveniently expressed as = ^xy. 

(b) Show that the net focal length of a pair of identical quadrupoles, one focusing 
and the next defocusing, and separated by a distance d, is 



OmBiL'^d' 

It is most practical to use transfer matrices which act on vectors (x x'), and first 
derive the transfer matrix for a drift space of length d and the focusing or defo- 
cusing quadrupole. The three matrices multiplied together gives the net transfer 
matrix, which contains the focal length. Assume the thin lens approximation. 

(c) Prove Equation IT?)1 

(d) Prove Equation IT9l in the limit that a « R and u « D. 

3. Suppose you are an experimenter on the MiniBooNE team. Following the example 
of |84j, you are about to place a temporary beam dump (see Figure EH) in the decay 
tunnel to stop a number of mesons and thereby test your beam flux prediction. By 
what fraction should the flux change? The z/g flux? Assume that half of z/g's come 
from K decay and half from /i decay for this experiment. Important information: the 
MiniBooNE decay pipe is 50 m, the temporary dump is at 25 m, and the focusing 
peak is ~800 MeV. 
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4. This problem concerns the distribution of neutrino daughters from pi/ii' decay. Assume 
that the meson parent is travelhng with (3 ^1. 

(a) Assuming the decay neutrinos are isotropic in the CM frame {dP/dQ' = l/47r=constant), 
show that in the laboratory frame 

dP/dn = {dP/dn') X (dn'/dn) = - , 

where 7 = (1 — /?^)~^/^ is the boost factor for the pion/kaon parent meson. In the 
small-angle limit tan^ << 1 this formula reduces to Equation El 

(b) Use your result above to show that the energy distribution in the lab of neutrinos 
from pion or kaon decays is flat in neutrino energy, i. e. 

dP/dEy = {dP/dn){dQ/dE^) = constant 

What is the constant? 
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